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L  REGIONAL  DISCRIMINATION  STUDIES 


1.1  INTRODUCTION 

Contract  #F19628-95-C-0184  for  "Regional  Discrimination  Studies"  was 
awarded  to  SMU  with  a  start  date  of  30  September  1995  and  ending  date  of  29 
September  1997.  Research  was  contingent  upon  the  use  of  data  from  the 
TXAR  array  that  was  installed  by  SMU  in  August  1993  under  Contract 
#F19628-93-C-0057.  Research  was  undertaken  to  calibrate  the  array,  and 
establish- a- "Ground  Truth  Data  Base."  It  soon  became  evident  that  acoustic 
data  was  needed  to  construct  such  a  data  base,  and  an  a  four-element 
infrasonic  array  was  installed  at  TXAR  in  February  1996. 

1.2  REPORTS  AND  OTHER  DELIVERABLES 

1.2.1  Reports 

Research  has  been  documented  in  Scientific  Reports  Nos.  1,  2,  and  3.  They 
consist  of  parts  as  described  below  plus  an  Acknowledgment  part  containing 
information  on  previous  contracts  and  publications  as  is  presented  in  Part  2 
of  this  report. 

1.2.1. 1  Scientific  Report  No.  1 

PL-TR-96-2249,  dated  April  1996,  consisted  of  three  parts.  Part  1  by  Chris 
Hayward  described  TXAR  in  detail.  Part  2  by  Heana  Tibuleac  and  Eugene 
Herrin  described  the  calibration  study  using  144  events  recorded  at  TXAR  for 
which  USGS  mb  values  were  available.  Part  3  by  G.  G.  Sorrells  and  Eugene 
Herrin  described  the  construction  of  the  ground  truth  data  base  using  a 
variety  of  sources  including  mine  explosions,  normal  earthquakes, 
earthquake  swarms,  and  very  shallow  earthquakes  induced  by  hydrocarbon 
production.  The  three  parts  are  summarized  in  the  following  sections. 
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1.2.1. 1.1  Part  1.  TXAR:  The  Prototype  Alpha  Array  Of  The  New  Improved 
Global  Seismic  System 

SMU  began  research  on  mini-array  technology  in  1991  under  Contract  #  MDA 
972-89-C-0054,  and  was  awarded  Contract  #  F19628-93-C-0057  to  install  TEXESS 
on  5  April  1993.  The  array  was  installed  by  SMU  personnel  at  their  LTX  site 
the  week  of  22  August  1993,  and  the  first  event  was  a  local  recorded  on  31 
August.  With  an  aperture  of  4  km,  TEXESS  [now  TXAR]  consists  of  nine 
sensor  sites,  which  includes  a  central  three-component,  short-period 
seismometer  installation  in  a  vault  at  the  hub,  and  eight  vertical  short-period 
seismometer  installations  in  surrounding  boreholes.  In  addition  to  the  short- 
period  instrumentation  at  the  hub,  a  posthole  K-S54000  long-period 
seismometer,  owned  by  SMU,  was  installed  in  a  shallow  borehole.  The  term 
posthole  was  coined  because  the  K-S54000  is  slimmer  than  the  K-S36000  and 
is  installed  in  a  shallow  hole  without  a  hole  lock  in  such  a  manner  that  it  can 
be  leveled  manually  rather  than  remotely.  Four  infrasonic  stations  were 
subsequently  established  using  porous  hose  arrays  and  inexpensive,  low- 
frequency  acoustical  sensors  collocated  with  the  seismometers  in  the 
boreholes.  The  TXAR  layout  including  infrasonic  stations  are  shown  in 
Figure  1-1  of  Appendix  1. 

TXAR  was  designed  as  an  improvement  upon  regional  arrays  of  the 
NORESS-ARCESS-GERESS  type.  Advancements  in  the  TXAR  design 
included  the  following: 

1.  The  placement  of  seismometers,  acoustical  sensors,  and  electronics  in 
boreholes  to  greatly  reduce  construction  costs  for  piers  and  vaults, 

2.  The  use  of  solar  power  at  each  site  rather  than  a  central-power  source, 

3.  The  use  of  GPS  receivers  for  time  data  at  each  seismometer  site  to  replace 
central  timing  from  the  hub, 

4.  The  employment  of  radio  links  from  seismometer  sites  to  the  hub  to 
replace  cable  links  and  associated  construction  costs, 

5.  The  use  of  modular  equipment  to  facilitate  the  installation  and 
maintenance  of  the  array. 
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The  above  advancements  in  array  design  reduced  costs  by  over  an  order  of 
magnitude  when  compared  with  GERESS.  Cost  per  element  at  GERESS  was 
$341,280  whereas  cost  per  element  at  TXAR  was  $33,000.  With  improved 
digital  instrumentation  to  assure  undistorted  phase  information  and  time- 
domain  processing,  azimuthal  deviations  for  test  events  have  been  reduced 
significantly.  As  a  result,  TXAR  is  better  than  either  a  single,  three- 
component  station  or  larger  arrays.  After  35  years,  we’re  back  to  the  original 
Geneva-type  array,  or,  in  the  words  of  Yogi  Berra,  "it's  dijd  vu  all  over  again.” 

1 .2.1.1 .2  Part  2.  Calibration  Studies  At  TXAR 

Calibration  studies  at  TXAR  used  a  modified  version  of  the  correlation 
method  described  by  Cansi,  Plantet  and  Massinon  in  order  to  estimate 
azimuth  and  horizontal  phase  velocity  of  144  events  for  which  USGS  mb 
values  were  available.  Modifications  to  the  correlation  method  include  the 
Fourier  interpolation  of  the  data  by  a  factor  of  8  to  obtain  a  virtual  sample  rate 
of  320 /sec,  use  of  an  LI  norm  (least  absolute  deviation)  to  obtain  estimates  of 
the  azimuth  and  phase  velocity  and  a  moving  window  display  to  indicate 
those  portions  of  the  waveform  that  show  strongest  correlation  across  the 
array. 

Corrected  phase  velocities  normally  associated  with  Pn  (less  than  8.6  km/s) 
are  generally  seen  for  events  at  epicentral  distance  as  far  as  2000  km.  For 
greater  distances,  upper  mantle  refracted  first  arrivals  (P)  with  corrected  phase 
velocities  greater  than  8.6  km/s  are  generally  observed  for  epicentral  distances 
beyond  1600  km.  Phase  identification  is  essential  in  order  to  select  a  suitable 
magnitude  scale. 

Based  on  the  144  well  located  events  (USGS)  and  using  the  Denny,  Taylor  and 
Vergino  formula,  the  most  reliable  magnitude  estimates  are  as  follows: 

1.  For  horizontal  phase  velocity  less  than  8.6  km/ sec: 

mb(D)=log  A  +  2.4  (log  D)  -  3.95  +  C  with  C=+0.3 

2.  For  horizontal  phase  velocity  greater  than  8.6  km/ sec 

mb(D)=log  A  +  2.4  Gog  D)  -  3.95  +  C  with  C=-0.50 
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The  M-discontinuity  beneath  TXAR  was  determined  to  the  first  order  to 
strike  along  an  azimuth  of  111  degrees  (NW-SE)  and  dip  10  degrees  to  the 
northeast.  This  result  is  consistent  with  the  tectonic  setting  for  the  area. 

1.2. 1.1. 3  Part  3.  Ground  Truth  Database 

The  preliminary  results  of  this  study  also  demonstrate  that  the  calculation  of 
a  running  estimate  of  the  normalized  correlation  coefficient  between  the 
outputs  of  collocated  seismic  and  infrasonic  sensors  provides  a  simple  but 
effective  method  for  the  detection  of  weak  short  period  acoustic  waves.  The 
success  of  this  approach  at  TXAR  is  attributed  to  the  fact  that  while  short 
period  acoustic  and  seismo-acoustic  signals  share  a  common  waveform,  short 
period  seismic  and  infrasonic  noise  are,  for  all  practical  purposes,  statistically 
uncorrelated  under  both  calm  and  windy  atmospheric  conditions. 
Statistically  uncorrelated  short  period  seismic  and  infrasonic  noise  is  expected 
to  be  a  property  of  all  "hard  rock"  geologic  environments  including  those 
which  are  likely  to  be  the  sites  of  future  IMS  installations.  However,  while 
acoustic  and  seismo-acoustic  signals  are  always  linearly  related,  they  will 
share  a  common  waveform  only  in  those  environments  where  the  seismic 
velocities  of  the  formation  containing  the  seismic  observation  point  are 
uniform  over  a  depth  range  that  is  large  compared  to  the  wavelengths  of  the 
acoustic  signals.  In  other  environments,  it  will  be  necessary  to  determine  the 
frequency  response  characteristics  of  the  seismo-acoustic  transfer  function  in 
order  to  optimize  the  performance  of  the  correlator  code.  In  this  regard,  it  has 
been  shown  that  the  existing  model  for  the  prediction  of  seismo-acoustic 
transfer  functions  suffered  from  serious  defects  when  it  was  tested  at  TXAR. 
Until  the  discrepancies  between  the  observed  and  predicted  values  of  the 
seismo-acoustic  scaling  factors  at  TXAR  can  be  explained,  experimental 
measurements  should  be  used  to  determine  the  frequency  response 
characteristics  of  local  seismo-acoustic  transfer  functions. 

1.2. 1.2  Scientific  Report  No.  2 
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FL-TR-96-2314,  dated  October  1996,  consisted  of  three  parts  and  an  app_udix. 
Part  1  by  Eugene  Herrin  and  G.  G.  Sorrells1  described  seismo-acoustic  studies 
at  TXAR.  Part  2  by  Jessie  L.  Bonner  desc  jed  a  study  of  surface  waves  recorded 
at  TXAR  from  quarry  explosions  near  Waco,  Texas.  Part  3  by  G.  G.  Sorrells  and 
Eugene  Herrin  described  the  use  of  acoustic  and  seismo-acoustic  data  to 
identify  vented  explosive  seismic  sources.  The  appendix  by  Chris  Hayward 
described  the  acoustical  sensors  installed  at  TXAR. 

1.2.1 .2.  1  Part  1.  Seismo- Acoustic  Studies  At  XAR 

The  use  of  seismic,  acoustic  and  seismo-acoustic  detections  has  been  clearly 
established  as  a  means  of  positively  identifying  surface  and  vented 
underground  explosions.  These  identified  events  then  become  part  of  a 
ground-truth  data  base.  This  procedure  can  be  used  at  collocated  seismic  and 
acoustic  arrays  that  detect  explosions  of  mb  2.5  or  larger  from  regional  sources. 

Detection  of  acoustic  signals  allows  the  identification  of  dusters  of  similar 
seismic  signals  from  explosions.  Under  light  wind  conditions  the  acoustic 
detection  threshold  is  a  few  tenths  of  a  microbar.  Under  moderate  wind 
conditions  the  seismo-acoustic  detection  threshold  is  about  a  microbar. 
Acoustic  signals  from  explosions  in  Coahuila,  Mexico,  with  seismic 
magnitudes  above  about  2.5  are  detectable  during  months  with  a  favorable 
zonal  stratospheric  wind  direction.  We  expect  to  detect  acoustic  signals  from 
explosions  in  New  Mexico  and  Arizona  this  winter  when  the  zonal 
stratospheric  winds  are  from  the  west. 

1.2.1.2.  2  Part  2.  Inversion  Of  Surface  Waves  For  Shallow  Velodty  Structure  In 
The  Fort  Worth  Basin 

The  Love  and  Rayleigh  wave  inversions  were  treated  as  separate  problems  to 
determine  if  the  upper  crust  in  the  Fort  Worth  basin  is  transversely 
anisotropic.  The  initial  results  show  two  different  velodty  structures  are 
obtained  through  the  separate  inversion  of  the  Rayleigh  (P-SV)  and  Love 


1  With  contributions  from  Jessie  L.  Bonner,  Valeriu  Burlacu,  Nancy  Cunningham,  Paul  Golden, 
Chris  Hayward,  Jack  Swanson,  Karl  Thomason,  and  Ileana  Tibuleac. 
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(SH)  waves  that  have  traveled  along  the  same  path.  This  is  not  unexpected 
considering  the  highly-  fractured  nature  of  the  rocks  in  the  upper  crust  of 
central  Texas  and  the  total  percentage  of  clays  in  the  stratigraphic  column. 
But  the  magnitude  of  this  phenomenon  is  unexpected  and  raises  doubts 
about  the  validity  of  the  results.  Further  research  will  work  to  alleviate  these 
doubts,  as  methods  for  inverting  Love  and  Rayleigh  waves  simultaneously 
using  code  that  considers  transversely  anisotropic  media  are  employed. 
Further  geophysical  data  from  the  basin  will  be  sought  to  constrain  the 
models  more  efficiently,  and  as  a  means  of  testing  the  final  answer. 

I.2.I.2.  3  Part  3.  A  Preliminary  Investigation  Of  The  Use  Of  Acoustic  And 
Seismo-Acoustic  Observations  To  Identify  Vented  Explosive  Seismic  Sources 

The  results  presented  above  demonstrate  that  short  period  acoustic  signals 
with  pressure  amplitudes  of  a  few  [ibars  or  less  will  generate  seismo-acoustic 
signals  that  are  detectable  at  a  seismically  "quiet"  site.  This  observation 
implies  that  they  will  also  be  detectable  at  future  IMS  sites  which  are  most 
likely  to  be  located  in  quiet  seismic  environments.  While  the  short  period 
acoustic  SNR  will  be  higher  at  the  outputs  of  the  infrasonic  monitoring 
system  during  periods  of  calm  to  light  surface  winds,  it  may  be  inferred  from 
the  results  that  the  seismo-acoustic  SNR  will  be  higher  during  periods  of 
moderate  to  high  surface  winds.  The  results  shown  also  indicate  that  seismic 
beamfonning  may  significantly  enhance  seismo-acoustic  SNR's  despite  the 
fact  that  the  seismic  array  is  not  optimally  configured  for  short  period  acoustic 
observations.  This  result  is  attributed  to  the  fact  that  unlike  seismic  signals, 
seismo-acoustic  signals  are  not  significantly  affected  by  multi-pathing  and 
remain  spatially  stable  over  the  characteristic  dimensions  of  a  typical  IMS 
seismic  array.  Thus,  it  may  be  inferred  that  the  seismic  monitoring  systems 
installed  at  future  IMS  sites  may  be  used  as  a  backup  infrasonic  system  during 
those  periods  when  local  atmospheric  turbulence  impairs  the  detection 
capability  of  the  primary  acoustic  monitoring  system. 

The  preliminary  results  of  this  study  also  demonstrate  that  the  calculation  of 
a  running  estimate  of  the  normalized  correlation  coefficient  between  the 
outputs  of  collocated  seismic  and  infrasonic  sensors  provides  a  simple  but 
effective  method  for  the  detection  of  weak-  short  period  acoustic  waves.  The 
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success  of  this  approach  at  TXAR  is  attributed  to  the  fact  that  while  short 
period  acoustic  and  seismo-acoustic  signals  share  a  common  waveform,  short 
period  seismic  and  infrasonic  noise  are,  for  all  practical  purposes,  statistically 
uncorrelated  under  both  calm  and  windy  atmospheric  conditions. 
Statistically  uncorrelated  short  period  seismic  and  infrasonic  noise  is  expected 
to  be  a  property  of  all  "hard  rock”  geologic  environments  including  those 
which  are  likely  to  be  the  sites  of  future  IMS  installations.  However,  while 
acoustic  and  seismo-acoustic  signals  are  always  linearly  related,  they  will 
share  a  common  waveform  only  in  those  environments  where  the  seismic 
velocities  of  the  formation  containing  the  seismic  observation  point  are 
uniform  over  a  depth  range  that  is  large  compared  to  the  wavelengths  of  the 
acoustic  signals.  In  other  environments,  it  will  be  necessary  to  determine  the 
frequency  response  characteristics  of  the  seismo-acoustic  transfer  function  in 
order  to  optimize  the  performance  of  the  correlator  code.  In  this  regard,  it  has 
been  shown  that  the  existing  model  for  the  prediction  of  seismo-acoustic 
transfer  functions  suffered  from  serious  defects  when  it  was  tested  at  TXAR. 
Until  the  discrepancies  between  the  observed  and  predicted  values  of  the 
seismo-acoustic  scaling  factors  at  TXAR  can  be  explained,  experimental 
measurements  should  be  used  to  determine  the  frequency  response 
characteristics  of  local  seismo-acoustic  transfer  functions. 

Finally,  the  results  of  this  study  infer  that  the  detection  of  a  short  period 
acoustic  signal  and  its  association  with  an  antecedent  near  regional  seismic 
event  is  a  promising  method  for  the  identification  of  vented  explosive 
sources.  If  further  research  confirms  this  premise  then  it  follows  that  the  data 
that  will  be  provided  by  the  collocated  seismic  and  infrasonic  monitoring 
systems  at  future  IMS  sites  may  be  used  to  significant  advantage  to  eliminate 
the  seismic  signals  generated  by  vented  explosions  from  further  consideration 
as  possible  indicators  of  a  clandestine  underground  nuclear  explosion. 

1.2.1 .3  Scientific  Report  No.  3 

PL-TR-97-2064,  dated  March  1997  consisted  of  two  parts.  Part  1  by  G.  G. 
Sorrells,  Eugene  Herrin,  and  Jessie  L.  Bonner2  described  the  construction  of  a 


2  With  contributions  from  Jack  Swanson,  Sarah  Deering,  and  Angela  Maddox. 
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regional  ground  truth  data  base  using  seismic  and  infrasonic  data  recorded  at 
TXAR.  Part  2  by  Jessie  L.  Bonner,  Sarah  Deering,  Tao  Liu,  Jack  Swanson,  and 
Heana  Tibuleac  described  seismic  and  infrasonic  data  observations  at  TXAR 
recorded  during  the  second  Seismo- Acoustic  Ground  Truth  (SAGII)  tests. 


1.2.1. 3.1  Part  1.  Construction  Of  Regional  Ground  Truth  Data  Bases  Using 
Seismic  And  Infrasound  Data 

The  results  of  the  TXAR  seismo-acoustic  experiment  illustrate  that 
infrasound  signals  generated  by  commercial  explosions  are  detectable  in  the 
0.5-5.0  hertz  bandwidth  and  can  be  observed  at  distances  up  at  least  680 
kilometers  from  the  source.  Furthermore,  it  may  be  inferred  from  these 
results  that  the  sources  of  the  majority  of  commercial  explosive  seismic 
events  will  be  identifiable  by  the  detection  of  associated  infrasound  signals 
when: 

•  the  magnitudes  of  the  seismic  events  are  greater  than  about  m^  =2.5 

•  their  source  to  receiver  distances  exceed  about  175-200  kilometers; 

•  the  events  occur  during  a  time  period  and  at  locations  such  that  the 
reception  point  is  located  stratospherically  "downwind"  from  the  sources. 

It  follows  then  that  application  of  the  seismo-acoustic  source  identification 
method  to  data  acquired  over  a  complete  stratospheric  zonal  wind  cycle  at  a 
particular  site  will  provide  the  information  required  to  construct  a  ground 
truth  data  base.  This  data  base  will  accurately  identify  the  sources  of  the 
seismic  events  that  occur  at  regional  distances  to  the  site  whose  magnitudes 
are  found  in  a  range  that  is  of  significance  to  the  verification  of  the  CTBT. 

1.2.1.3.2  Part  2.  Seismic  And  Infrasound  Data  Observations  At  TXAR 

The  second  seismic  and  acoustic  ground  truth  test  (SAGII)  using  the  TXAR 
array  was  initiated  on  1  February  1997.  A  goal  of  this  test  is  to  determine  the 
percentage  of  seismic  events  recorded  at  TXAR  between  November  1- 
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December  15,  1996  that  had  associated  acoustic  signals.  During  this  time 
period,  stratospheric  zonal  winds  blew  from  west  to  east,  thus  enhancing  the 
opportunity  to  record  infrasound  signals  from  the  west  of  TXAR.  This 
temporal  factor  also  decreased  the  chance  of  recording  infrasound  signals 
originating  east  of  the  array,  an  inference  confirmed  by  the  preliminary 
results  of  SAGTH  To  date,  infrasound  signals  have  only  been  associated  with 
regional  events  from  the  west  of  the  array  for  a  time  period  between  1  and  7 
November  1996. 

1.2.2  Other  Deliverables 

Software  developed  for  the  contract  will  be  delivered  in  accordance  with  the 
Contract  Data  Requirements  List  (CDRL). 
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NORSAR,  and  Provide  the  Interface  to  the  Intelligent  Array  System  at 
NORSAR. 

Quarterly  R&D  Status  Report,  4  Jan.  to  15  Apr.  1989. 

Quarterly  R&D  Status  Report,  15  Apr.  to  15  July  1989. 

Quarterly  R&D  Status  Report,  15  July  to  15  Oct.  1889. 

Quarterly  R&D  Status  Report,  15  Oct.  1889  to  15  Jan.  1990. 

Quarterly  R&D  Status  Report,  15  Jan.  to  15  Apr.  1990. 

Quarterly  R&D  Status  Report,  15  Apr.  to  15  July  1990. 
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APPENDIX  1.  SEISMO-ACOUSTIC  STUDIES  AT  TXAR1 

Eugene  Herrin  and  G.  G.  Sorrells 
with  contributions  from 

Jessie  Bonner,  Nancy  Cunningham,  Sarah  Deering,  Paul  Golden,  Chris  Hayward,  Tao 
Liu,  Angela  Maddox,  Jack  Swanson,  Karl  Thomason,  and  Deana  Tibuleac 


1.1  INTRODUCTION 

TXAR  is  located  in  far  West  Texas  only  a  few  miles  north  of  the  border  with  Mexico. 
The  seismic  sources  for  regional  events  recorded  at  TXAR  are  scattered  mainly 
throughout  northern  Mexico,  Texas,  New  Mexico,  and  Arizona.  Early  analysis  indicated 
that  at  least  15%  of  these  events  were  commercial  explosions.  Whereas  the  mine 
operators  in  the  U.S.  have  been  helpful  in  identifying  their  explosions,  we  have  not  been 
able  to  obtain  similar  information  from  mines  in  Mexico.  We  established  a  four  element 
infrasonic  array  within  TXAR  so  that  infrasound  signals  from  surface  explosions  could  be 
associated  with  seismic  data  in  order  to  identify  mining  explosions,  particularly  in 
Mexico  (see  Figure  1)..  This  method  has  proved  to  be  very  useful  during  periods  when 
stratospheric  winds  are  in  the  source-receiver  direction.  Identification  of  a  variety  of 
seismic  sources  is  possible  using  the  seismo-acoustic  data.  These  include  mining  and 
construction  explosions,  local  explosions,  shuttle  quakes,  and  bolide  quakes.  This  poster 
illustrates  techniques  used  to  construct  a  ground-truth  database  at  TXAR  for  the 
Southwestern  United  States. 


1  Herrin,  Eugene,  1997,  Ground  truth  for  regional  seismic  events  recorded  at  TXAR:  Abstracts  [poster]  of 
the  19th  Annual  Seismic  Research  Symposium  on  Monitoring  a  Comprehensive  Test  Ban  Treaty,  30  August 
1997,  p.  37. 
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Figure  1-1.  TXAR  location  map.  Acoustic  sensors  are  located  at  TXOO  TX02  TX03 
and  TX09.  ’  ’ 
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Figure  1-2.  Locations  for  events  recorded  at  TXAR  between  July  15,  1996  and 
August  10,  1996.  The  majority  of  the  epicenters  occur  in  four  clusters  within  or 
near  known  mining  complexes.  The  cluster  found  east-southeast  of  TXAR  at  a 
range  of  300-320  km  is  approximately  coincident  with  the  MI  CARE  (Minera  de 
Carbonifera  Rio  Escondito)  coal  mining  district  located  in  Coahuila,  Mexico  near 
the  border  with  the  United  States.  The  cluster  found  due  South  of  TXAR  at  a 
range  of  about  140-160  km  occurs  near  the  Minas  de  Hercules  iron  mines.  The  two 
clusters  found  about  550-700  km  northwest  of  TXAR  coincide  with  the  Morenci- 
Silver  City  Porphyry  copper  mining  district  in  southeastern  Arizona  and 
southwestern  New  Mexico. 
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Local  Explosion--  February  11, 1997 


Figure  1-4.  Seismic  event  on  February  11,  1997  originating  within  10  km  of  the 
TXAR  array.  For  several  years,  similar  events  were  thought  to  be  microseismic 
activity  along  nearby  faults.  Implementation  of  the  TXAR  infrasound  array  in 
March,  1996,  allowed  us  to  determine  the  origin  of  these  events,  with  duration 
magnitudes  <  -0.5,  to  be  explosions  southwest  of  Terlingua,  Texas.  The  amplitudes 
are  plotted  at  a  uniform  scale  to  show  the  attenuation  of  the  signal  across  the  array. 
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Figure  1-5.  Seismograms  (A)  and  microbarograms  (B)  from  an  explosion  in  the 
MI  CARE  mining  district  of  northern  Mexico. 
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Figure  1-6.  Balanced  zonal  winds  at  30°N  latitude.  Wind  from  the  west  (positive)  is 
indicated  by  light  colors,  while  winds  from  the  east  (negative)  are  shown  by  dark 
colors. 
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Figure  1-7.  A)  Seismograms  generated  by  an  explosion  at  the  Tyrone  open-pit 
copper  mine  in  Silver  City,  New  Mexico.  B)  Microbarograms  of  the  infrasound 
signal  associated  with  reflection  from  the  upper  atmosphere.  The  group  velocity  of 
this  signal  is  0.3  km/s  referenced  to  the  origin  time  and  location  of  the  explosion. 
The  microbarograms  have  been  time  aligned  for  an  horizontal  phase  velocity  of  0.34 
km/s  and  a  back-azimuth  of  306. 
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Figure  1*8.  Vertical-component  seismogram  from  a  explosion  (with  no  delays 
between  rows  or  holes)  at  the  Phelps-Dodge  mine  at  Tyrone,  New  Mexico.  The  inset 
shows  three  distinct  P  wave  arrivals.  The  second  of  these  arrivals  is  thought  to  be 
the  seldom  observed  phase,  PnPn.  This  phase  has  not  been  observed  on  records 
from  ripple-fired  blasts  originating  at  the  Tyrone  mine,  suggesting  that  the  PnPn 
arrival  is  obscured  by  the  coda  of  Pn  for  these  explosions.  Arrival  times  for  the 
phase  are  in  agreement  with  those  predicted  by  the  IASPI  model.  These  explosions 
are  more  likely  to  produce  infrasound  signals  seen  at  regional  distances  than  ripple- 
fired  explosions  from  the  same  mine. 
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Tyrone  AZ  explosion  1 1  October  1996  546.0  km  0.75-5  Hz  shifted:  325  deg  0.35  km/s 
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Figure  1-9.  Results  of  semblance  processing  on  seismo-acoustic  data  from  the 
October  11, 1996  Tyrone  explosion. 
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1.2  SHUTTLE  QUAKE  SUMMARY 


The  location  of  the  TXAR  array  allows  us  to  record  the  N  waves  produced  as  the  space 
shuttle  makes  its  landing  approach  to  Kennedy  Spaceflight  Center.  While  the  N  waves 
recorded  from  the  Shuttle  are  spectacular  (-100  microbars),  the  unusual  components  of 
these  signals  were  recorded  on  TXAR’s  seismic  channels.  Short-period  seismic  waves 
are  seen  prior  to  the  N  wave  arrival  that  could  easily  be  mistaken  for  surface  explosions. 
Phases  recognized  in  these  signals  include  Pg,  Sg,  and  Rg.  FK  processing  of  the  phases 
showed  appropriate  phase  velocities  but  different  back-azimuths  for  each  phase.  We  have 
termed  these  small  seismic  events  “shuttle  quakes”,  and  the  results  of  our  initial  study  are 
presented  in  this  section  of  the  poster  presentation. 
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STS82 


Time  (seconds)  From  08:37:06  UTC  on  February  21, 1997 


Figure  1-10.  Seismograms  and  barograms  from  the  STS-82  shuttle  mission. 
Distinct  phases  include  Pg  and  short-period  Rayleigh  (Rg).  The  infrasound  sensors 
of  the  TXAR  seismo-acoustic  array  recorded  the  supersonic  pressure  wave  (N 
wave). 
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STS-82  for  1 87.5  degrees  azimuth  -  level  2,  corrected  for  elevation 


Figure  1-11.  Results  of  processing  the  N  wave  using  wavelet  decomposition. 
Horizontal  phase  velocity  and  back-azimuth  were  determined  to  be  0.72  km/s  and 
188.8,  respectively.  Frequency-wavenumber  (F-K)  analysis  was  used  to  determine 
the  apparent  phase  velocities  and  back-azimuths  for  the  seismic  phases. 
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STS75 


Time  (seconds)  From  13:43:02  UTC  on  March  9, 1996 


Time  (seconds)  From  13:43:02  UTC  on  March  9,  1996 


Figure  1-13.  Shuttle  quake  from  the  STS-75  mission  which  passed  approximately 
150  km  north  of  the  TXAR  array. 
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Seismic 


STS70 


Figure  1*14.  Shuttle  quake  from  the  STS*70  mission  which  passed  approximately  50 
km  south  of  the  TXAR  array. 
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Time  (seconds)  From  11:22:34  UTC  on  Sept.  18, 1995 


Figure  1-15.  Shuttle  quake  and  N  wave  from  the  STS-69  mission  which  passed 
approximately  150  km  north  of  the  TXAR  array. 
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STS80 


Time  (seconds)  From  11:44:49  UTC  on  Dec.  7, 1996 


Figure  1-16.  Infrasound  from  the  STS-80  mission  which  passed  far  enough  from 
TXAR  that  no  N  wave  or  shuttle  quake  was  observed. 
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STS-82  SHUTTLE  QUAKE 

February  21, 1997  Landing  Time:  8:32:26  UTC 
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Figure  1-17.  Initial  results  of  the  “Bow  Wave”  theory  using  phase  velocity  and  back- 
azimuth  data  from  STS-82.  Results  are  promising,  and  additional  shuttle  quake 
data  are  now  being  used  to  further  test  the  theory. 
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Figure  1-18.  Wavelet  transform  of  the  N  wave  from  STS-82.  The  arrival  of  the  N 
wave  on  the  seismic  channels  is  seen  as  a  sharp  impulse  on  Levels  1-4  of  the 
transform. 
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STS -82  waveforms  at  TX10 


Discrete  Wavelet  Transform,  level  5,  using  Haar  wavelet,  absolute  coefficients  at  TX10 


Discrete'Wavefet  Transform,  level  5,  using  Haar  wavelet,  absolute  coefficients  at  TX08 


Sample  points 


Figure  1-19.  Wavelet  decomposition  of  the  N  wave  at  TX01.  The  arrival  time  for  the 
phase  is  established  to  better  than  one  sample  point. 
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1.3  BOLIDE  SUMMARY 


The  Geminid  meteor  shower  is  active  between  December  6  and  December  19,  reaching 
peak  activity  on  December  13.  Unlike  most  meteor  showers  which  are  by-products  of 
comets,  the  Geminid  orbit  does  not  match  any  known  orbit  for  a  comet.  However,  the 
Geminid  orbit  does  match  the  Earth-crossing  asteroid  3200  Phaethon,  a  four-mile  wide 
object  discovered  in  1983.  Thus  the  Geminids  are  the  first  meteor  shower  to  be  positivelv 
linked  to  asteroids.  Rates  during  the  peak  observation  time  are  80  per  hour.  SMU 
personnel  at  TXAR  between  December  9  and  December  13  reported  sighting  many 
meteors.  On  December  12,  1996,  we  recorded  over  120  acoustic  arrivals  from  the 
northwest  of  TXAR.  Based  upon  the  location  of  the  radiant  of  the  Geminids,  we  believe 
these  signals  were  the  result  of  the  meteors  entering  the  upper  atmosphere.  This  final 
section  of  the  poster  presentation  looks  at  suspected  bolide  signals  at  TXAR. 
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Figure  1-20.  Arrival  times  for  127  infrasound  signals  recorded  at  TXAR  on 
December  12, 1996  plotted  against  back-azimuth.  These  signals  are  suspected  to  be 
from  meteors  of  the  Geminid  shower  entering  and  causing  acoustic  disturbances  in 
the  upper  atmosphere  northwest  of  TXAR  (290-340  degrees).  The  maximum 
amplitude  of  this  signal  is  approximately  40  mbars,  one  of  the  largest  infrasound 
events  recorded  at  TXAR.  The  gap  in  the  data  between  2:00  and  14:00  GMT  could 
be  related  to  the  zonal  wind  effect  discussed  earlier.  Half  of  the  day’s  events  should 
be  in  the  atmosphere  northeast  of  TXAR,  and  thus  would  be  attenuated  by  strong 
winds  in  the  upper  atmosphere. 
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Geminkd  event  12  December  1996  0.75-5  Hr  shifted:  337  deg  0.35  km/s 
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Figure  1-21.  Semblance  between  traces  for  a  Geminid  acoustic  event  on  December 
12,1996. 


APPENDIX  2.  OBSERVATIONS  OF  SHUTTLE  QUAKES  AT  TXAR 


Jessie  L.  Bonner 
Eugene  T.  Herrin 
Gordon  G.  Sorrells 
Jack  G.  Swanson 
Heana  M.  Tibuleac 

Introduction 

People  in  the  remote  ghost  town  of  Terlingua,  Texas,  were  awakened  early  on 
February  21, '1997  by  a  “large  explosion”.  SMU  Geophysics  operates  a  ten-element 
seismic  and  four-element  acoustic  array  in  this  area  and  was  contacted  by  a  concerned 
citizen  and  asked  to  help  identify  the  cause  of  the  large  acoustic  event.  After  looking  at 
the  classic  “N  wave”  on  the  acoustic  channels,  we  concluded  the  residents  had  heard  a 
sonic  boom.  The  name  “N  wave”  comes  from  the  characteristic  sharp  pressure  increase 
followed  by  a  linear  decrease  in  pressure  until  ambient  pressure  is  restored  by  a  sharp 
pressure  rise  (see  Figure  1).  After  checking  the  NASA  web-pages,  we  knew  the  sonic 
boom  was  produced  from  the  space  shuttle,  which  had  landed  at  Kennedy  Space  Center 
some  15  minutes  after  the  explosion  had  been  heard.  By  looking  at  the  landing 
groundtrack,  we  concluded  that  shuttle  missions  STS-82  had  passed  less  than  55  km 
south  of  the  TXAR  array  on  its  landing  approach  to  the  Kennedy  Spaceflight  Center. 

While  the  N  waves  from  STS-82  were  spectacular  (-100  microbars),  the  truly 
remarkable  parts  of  these  signals  were  recorded  on  the  seismic  channels  of  the  TXAR 
array.  By  looking  at  Figure  1,  one  might  conclude  that  the  sonic  wave  arrived  and  was 
recorded  at  TXAR  during  another  seismic  event.  The  presence  of  short-period, 
fundamental  mode  Rayleigh  waves,  Rg,  in  this  seismic  signal  suggests  a  very  shallow 
source.  However,  initial  FK  processing  of  the  Pg  and  Rg  waves  showed  appropriate 
phase  velocities  but  different  back-azimuths  (Table  I)  for  each  phase. 

Gordon  Sorrells  and  Eugene  Herrin  then  proposed  that  these  seismic  signals  were 
“bow  waves”  created  as  the  supersonic  N  wave  sweeps  across  the  surface  of  the  earth. 
Thus,  work  started  on  the  theory,  and  results  produced  thus  far  have  been  promising. 
Predicted  back-azimuths  for  different  phases  (Pg  and  Rg)  are  within  5°  of  the  true  back- 
azimuths.  We  have  termed  these  small  seismic  events  created  by  the  passage  of  the 
shuttle  as  “shuttle  quakes”. 

This  paper  represents  a  compilation  of  the  data  recorded  thus  far  at  TXAR  from 
passage  of  different  shuttle  missions.  Appendix  3  shows  work  completed  on  the  theory. 
So  far,  seven  different  shuttle  missions  have  had  landing  groundtracks  close  enough  to 
TXAR  for  recording  of  both  shuttle  quakes  and  N  waves,  or  infrasound  signals.  For  each 
mission,  we  show  plots  of  the  data  as  well  as  tables  showing  parameters  for  each  mission 
(Figs.  1-18  and  Tables  I-IX).  These  observations  will  be  combined  with  the  “bow  quake” 
theory,  as  well  as  the  wavelet  analysis  results,  and  submitted  to  Seismological  Research 
Letters  within  the  upcoming  months. 
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Processing  Seismo-Acoustic  Data 


Horizontal  phase  velocities  and  back-azimuths  for  each  of  the  phases  identified  on 
the  “shuttle  quakes”  were  determined  by  frequency-wavenumber  (F-K)  analysis.  This 
technique  transforms  the  array  data  from  a  time-space  domain  into  a  frequency- 
wavenumber  domain.  We  applied  a  2-second  window  around  each  phase  (P,  S,  and  Rg 
arrivals)  and  used  the  F-K  technique  to  determine  back-azimuth  and  phase  velocity.  For 
body  wave  arrivals,  the  data  was  filtered  with  a  3rd  order  bandpass  filter  between  0.6  and 
4.5  Hz  to  optimize  results.  For  surface  waves,  this  was  modified  to  a  3rd  order  bandpass 
filter  between  0.1  and  2  Hz.  Results  from  previous  studies  that  used  regional  seismic 
events  have  shown  the  F-K  technique  to  be  very  reliable  at  TXAR  with  standard 
deviations  'of  approximately  ±0.2  km/sec  and  ±5  degrees  for  phase  velocity  and  back- 
azimuth,  respectively.  The  results  of  the  F-K  processing  are  shown  in  Tables  I,  II  El,  V, 
and  VI.  Notice  how  the  back-azimuths  change  from  phase  to  phase.  The  same  technique 
was  used  to  determine  back-azimuth  for  the  infrasound  arrivals  from  STS-78  and  STS-80 
(Table  VIII  and  IX),  both  missions  being  at  great  distances  from  TXAR. 

A  different  technique  was  used  to  determine  the  back-azimuth  and  phase  velocity 
for  the  supersonic  “N  waves”.  Since  the  TXAR  infrasound  array  has  only  been  in 
operation  since  March  1996,  we  could  not  determine  azimuths  for  the  N  waves  using  the 
infrasound  channels  for  missions  prior  to  STS-78.  Instead,  we  chose  to  use  the  seismic 
response  to  the  N  wave  as  recorded  on  the  vertical-component  seismometers.  The 
seismic  response  in  dominated  by  two  pressure  increases  at  the  beginning  and  end  of  the 
boom  which  produce  downward  spikes  on  the  seismograms.  However,  determining 
where  these  pressure  increases  began  proved  problematic  as  the  “shuttle  quakes”  masked 
the  N  wave  arrival  for  missions  STS-82  and  STS-70. 

Thus,  we  decided  to  use  wavelet  analysis  for  the  problem.  The  Haar  wavelet, 
which  is  discontinuous  and  resembles  a  step  function,  was  considered  optimal  for 
detection  of  sharp  discontinuities  in  a  signal  otherwise  masked  by  the  simultaneous 
seismic  arrivals.  The  arrival  times  of  the  N  waves  on  the  seismic  channeis  were  then 
obtained  with  a  precision  of  a  sample  point  by  discrete  wavelet  decomposition,  and  then 
used  to  calculate  the  horizontal  phase  velocity  and  back-azimuth.  With  one  exception, 
see  Table  IV,  the  second  level  of  detail  in  the  wavelet  decomposition  for  each  N  wave 
was  used  to  manually  find  the  arrival  times  at  each  station  in  the  TXAR  array.  Since  a  N 
wave  has  a  spectral  response  at  all  frequencies,  the  distribution  of  energy  in  all  the  other 
levels  were  considered  as  well.  Well  defined  N  waves  from  events  to  the  south  of  TXAR 
(STS-70  and  STS-82)  worked  best  in  the  decomposition  routine.  N-waves  from  shuttle 
passages  North  of  TXAR  (STS-69,  STS-78,  and  STS-92)  were  less  defined,  but  because 
of  refraction  and  reflection  of  infrasound  in  the  atmosphere,  more  than  one  arrival  was 
processed  (see  Tables  IV  and  VH). 
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Figure  2-1.  Seismograms  and  barograms  from  the  STS-82  shuttle  mission.  Distinct 
phases  include  Pg  and  short-period  Rayleigh  (Rg)  recorded  on  the  vertical 
component  of  the  seismometers  of  the  TXAR  array.  The  infrasound  sensors  of  the 
TXAR  seismo-acoustic  array  recorded  the  classic  “N  wave”  from  the  passage  of  the 
supersonic  pressure  wave.  This  wave  coupled  into  the  ground  and  can  be  seen  as  a 
packet  of  different  frequencies  on  the  seismic  channels  (labeled  as  N  on  the  upper 
plot). 


Time  (seconds)  From  08:37:06  UTC  on  February  21 , 1997 


STS-82  PARAMETERS 


Shuttle:  Discovery 

Mission  Dates:  February  11-21,  1997 

Landing:  KSC  at  08:32:26  UTC 

Mission:  Repair  mission  II  for  the  Hubble  space  telescope 

Landing  Groundtrack:  Approximately  55  km  south  of  TXAR  bearing  80° 


KSQ  2nd  LANDING  OPPORTUNITY 

Deorhrt  Bum:  am  CST 


Landing: 


9:27  am  CST 


Tiroes  do  not  correspond  to  STS-82  landing. 


Table  I.  Mission  information  and  results  of  processing  seismic 
STS-82. 


and  acoustic  signals  from 


1  Measured  at  station  TX06 

'  Based  upon  convergence  of  FK  results  to  a  stable  solution. 

3  Short-period  Rayleigh  waves  caused  by  coupling  of  the  N  wave  into  the  near-surface.  Since  this  phase  has 
a  velocity  greater  than  the  N  wave,  the  phase  runs  out  in  front  of  the  N  wave.  However,  determining  the 
back-azimuth  and  phase  velocity  is  problematic  due  to  the  lack  of  coherence  in  the  phase  across  themrav 
Obtained  using  a  wavelet  2  decomposition  of  the  seismo-acoustic  channels. 
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STS-82  for  187.5  degrees  azimuth  -  level  2.  corrected  for  elevation 


Distance  (km)  relative  to  TX06  in  the  propagation  direction 

Figure  2-2.  Results  of  processing  the  “N  wave”  for  STS-82  using  wavelet 
decomposition.  These  results  are  corrected  for  elevation  differences  of  the  stations 
at  TXAR. 


STS-82  for  187.1  degrees  azimuth  -  level  2 


Distance  (tan)  relative  to  TX06  in  the  propagation  direction 

Figure  3-3.  Results  of  processing  the  “N  wave”  for  STS-82  using  wavelet 
decomposition.  These  results  are  not  corrected  for  elevation  differences  of  the 
stations  at  TXAR. 


STS70 


Figure  2-4.  Seismograms  produced  from  the  shuttle  mission  STS-70.  We  checked 
our  archived  array  data  for  other  shuttle  events,  and  found  seven  separate  missions 
which  created  either  seismic  or  infrasound  that  was  recorded  at  TXAR.  STS-70  had 
a  similar  groundtrack  to  STS-82,  thus  back-azimuths  and  phase  velocities  (Table  II) 
were  similar.  In  fact,  spectral  comparisons  show  how  similar  the  Pg  and  Rg  signals 
are  for  STS-82  and  STS-70  (Figs.  7-8). 
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STS-70  PARAMETERS 


Shuttle:  Discovery 

Mission  Dates:  July  13-22,  1995 

Landing:  KSC  at  12:02:11  UTC 

Mission:  Deployment  of  the  7th  Tracking  Data  and  Relay  satellite 
Landing  Groundtrack:  .Approximately  55  km  south  of  TXAR  bearing  80° 


KSC  2ND  LANDING  OPPORTUNITY 

Deorbit  Bum:  am  CST 

Landing:  ^:27  am  CST 


Tiroes  do  not  correspond  to  STS-70  landing. 


Phase 

Arrival  Time5 

Back-azimuth 

Velocity 

Quality6 

11:47:40.73 

.63 : 

11:47:53.18 

-r-5r2.5'-  - . 

11:48:06.43 

!  ‘ 1  -  ’If  ~~  -LI.  Ti.-J 

2.4 

11:48:10 

Table  H  Mission  information  and  results  of  processing  seismic  and  acoustic  signals  from 
STS-70. 


3  Measured  at  station  TX06 
6  Based  upon  FK  convergence  to  a  stable  solution 

'  No  acoustic  channels,  thus  measurements  obtained  using  seismo-acoustics  and  a  level  2  wavelet 
decomposition. 
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Distance  (km)  relative  to  TX06  in  the  propagation  direction 


Figure  2-5.  Results  of  processing  the  “N  wave”  for  STS-70  using  wavelet 
decomposition.  These  results  are  corrected  for  elevation  differences  of  the  stations 
atTXAR. 

STS-70  for  185.8  degrees  azimutfi  -  level  2 


Figure  2-6.  Results  of  processing  the  “N  wave”  for  STS-70  using  wavelet 
decomposition.  These  results  are  not  corrected  for  elevation  differences  of  the 
stations  at  TXAR. 
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Figure  2-7.  TX02  seismograms  from  the  STS-82  and  STS-70  missions  show  how 
similar  the  two  shuttle  quakes  are.  The  increased  separation  in  the  phases  in  the 
lower  plot  suggests  that  STS-70  flew  further  south  of  TXAR  than  STS-82. 
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Frequency  (Hz) 


Figure  2-8.  Pg  and  Rg  spectra  for  STS-82  and  STS-70. 
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STS69 


Figure  2-9.  Shuttle  quake  from  the  STS-69  mission  which  passed  approximately  150 
km  north  of  the  TXAR  array. 
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STS-69  PARAMETERS 


Shuttle:  Endeavor 

Mission  Dates:  September  7-18,  1995 
Landing:  KSC  at  1 1:37:56  UTC 

Mission:  Experiments  with  the  Wake  Shield  Facility  and  deployment  of  the  Spartan  201 
astronomy  satellite 

Landing  Groundtrack:  Approximately  150  km  north  of  TXAR  bearing  95° 


KSG  1ST  LANDING  OPPORTUNITY 

Denrbit  Bum:  6:49  am  CS7 
Landing:  7:52  am  CST 


Times  do  oot  correspond  to  STS-69  landing. 


Phase 


Arrival  Time8 

Back-Azimuth 

Velocity 

11:23:3535 

WFiim 

11:23:45.53 

11:24:05.03 

i  vy  33 

.11:24:25.85 

1.,,  2.1  •  •; 

Table  Ed.  Mission  information  and  results  of  processing  seismic  signals  from  STS-69. 


8  Measured  at  TX06 
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Table  TV.  Wavelet  processing  of  inffasound  arrivals  from  STS-69. 


STS-69  for  348.9  degress  azimuth  -  first  arrival  level  2 


Figure  2-10.  Results  of  processing  the  first  infrasound  arrival  for  STS-69  using 
wavelet  decomposition. 


Figure  2-11.  Results  of  processing  the  second  infrasound  arrival  for  STS-69  using 
wavelet  decomposition. 


Figure  2-12.  Results  of  processing  the  third  infrasound  arrival  for  STS-69  using 
wavelet  decomposition. 
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Acoustic 


STS75 


Time  (seconds)  From  13:43:02  UTC  on  March  9, 1996 


Figure  2-13.  Shuttle  quakes  and  N  waves  from  STS-75. 
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STS-75  PARAMETERS 


Shuttle:  Columbia 

Mission  Dates:  February  22-March  9,  1996 
Landing:  KSC  at  13:58:38  UTC 

Mission:  Deployment  of  the  Tethered  Satellite  System  Reflight 
Landing  Groundtrack:  Approximately  150  km  north  of  TXAR  bearing  95° 


iphase  «^rrival  Time9 10 


13:44:06.95  S 

13:44:15.1 

13:44:32.7 

13:44:54.7 

13:45:19.6 

13:48:273 


Table  V.  Mission  information  and  results  of  processing  seismic  and  acoustic 
STS-75. 


signals  from 


9  Measured  at  TX06 

10  Only  two  acoustic  channels  available,  thus  measurements  obtained  using  level  2  wavelet  decomposition 
of  the  seismo-acoustics  channels. 
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12 


STS-75  for  349.9  degrees  azimuth  -  level  2.  corrected  for  elevation 


Figure  2-14.  Results  of  processing  the  infrasound  arrival  for  STS-75  using  wavelet 
decomposition.  The  results  are  corrected  for  elevation  differences  of  stations  at 
TXAR. 


STS-75  for  349.8  degrees  azimuth  -  level  2 


Figure  2-15.  Results  of  processing  the  infrasound  arrival  for  STS-75  using  wavelet 
decomposition.  The  results  are  not  corrected  for  elevation  differences  of  stations  at 


STS94 


Time  (seconds)  From  10:31:31.5  UTC  on  July  17, 1997 


Time  (seconds)  From  10:31:31.5  UTC  on  July  17, 1997 


Figure  2-16.  Shuttle  quakes  and  N  waves  from  STS-94. 
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STS-94  PARAMETERS 


Shuttle:  Columbia 

Mission  Dates:  July  1-17,  1997 

Landing:  KSC  at  10:46:34  UTC 

Mission:  Reflight  of  Microgravity  Science  Laboratory 

Landing  Groundtrack:.  Approximately  150  km  north  of  TXAR  bearing  80°. 


Arrival  Time11 


10:32:09.65  > 
10:32:133 
I  10:32:41.13 
10:33:03.9  J 


Rack-azimuth 


Table  VI.  Mission  information  and  results  of  processing  seismic  signals  from  STS-94. 


Arrival  #  Vh  Back-azimuth 


,:V~348.7* 

SgwTilgf  ’  v 

"V>-:  3453 •• 


Table  VH.  Wavelet  processing  of  infrasound  arrivals  from  STS-94. 


1 1  Measured  at  TX06 
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STS- 94  for  345.3  degrees  azimuth  -  third  arrival,  level  2,  corrected  for  elevation 


Figure  2-19.  Wavelet  processing  of  the  third  infrasound  arrival  from  STS-94. 
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STS78 


o 

o? 

3 

O 

a 

< 


Time  (seconds)  From  12:42:13  UTC  on  July  7, 1996 

Figure  2-20.  Infrasound  from  the  shuttle  mission  STS-78  recorded  at  TXAR.  The 
shuttle  passed  far  enough  from  TXAR  that  no  N  wave  or  shuttle  quake  was 
observed. 
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STS-78 


Shuttle:  Columbia 

Mission  Dates:  June  20-July  7, 1996 

Landing:  KSC  at  12:37:30  UTC 

Mission:  Experimenting  with  the  Life  and  Microgravity  Spacelab 
Note:  No  shuttle  quake 
Landing  Groundtrack: 


KSC  1ST  LANDING  OPPORTUNITY 

<>:IW  am  CDT 

% 


Times  do  not  correspond  to  STS-78  landing. 


Table  VHL  Mission  information  and  results  of  processing  acoustic  signals  from  STS-78. 
Five  different  acoustic  signals  were  processed. 


12  Measured  from  TXI03 
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Acoustic 


STS80 


Figure  2-21.  Infrasound  from  the  shuttle  mission  STS-80. 
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STS-80 


Shuttle:  Columbia 

Mission  Dates:  November  19-December  7,  1996 
Landing:  KSC  at  11:49:05  UTC 
Mission:  3rd  Flight  of. the  Wake  Shield  Facility 
Note:  No  shuttle  quake 
Landing  Groundtrack: 


j  -  Phase 

Arrival  Time13 

Azimuth 

Velocity  Quality  | 

;  41:46:57.6  >  ; 
^41:47:15.7  c;: 
'741:47:44.6 
.11:48:454 3 
:  11:50:00.9  * 

7i n 

-?mm- 

msi 

Table  DC.  Mission  information  and  results  of  processing  acoustic  signals  from  S  l  S-80. 
Five  separate  acoustic  signals  were  processed. 


13  Measured  from  TXI03 
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APPENDIX  3.  SEISMIC  "BOW"  WAVES 
G.  G.  Sorrells 

PART  I.  FORMULATION  OF  THE  PROBLEM 

Objects  such  as  meteorites  and  space  shuttle  craft  can  travel  at  speeds  in  excess  of  mach 
10.  Since  maximum  shock  overpressures  scale  as  the  square  of  the  mach  number,  a 
significant  shock  front  will  be  developed  and  sustained  by  the  movement  of  these  objects 
through  the  atmosphere.  The  intersection  of  the  shock  front  with  the  surface  of  the  earth 
will  generate  elastic  disturbances  which,  depending  upon  the  ground  speed  of  the  shock 
front  relative- to  the  local  horizontal  seismic  phase  velocities  will  be  seen  by  observers 
near  the  ground  path  of  the  object  as  either  propagating  seismic  waves,  or  as 
inhomogeneous  waves.  These  waves  are  referred  to  in  this  paper  as  seismic  "bow"  waves 
because  of  their  resemblance  to  the  fluid  waves  generated  by  the  movement  of  a  ship's 
bow  through  water.  An  approximate  theory  which  explains  the  generation  of  seismic 
"bow"  waves  by  a  hypersonic  shock  front  is  developed  in  the  following  paragraphs. 


Figure  3-1.  Trace  of  an  idealized  shock  front  in  the  X-Y  plane  at  time  t\  Position 
vectors  necessary  to  describe  the  seismic  response  to  the  shock  front  are  shown  as 
arrows  (See  text  for  additional  details) 
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Definitions 

The  surface  of  the  earth  is  approximated  by  the  X-Y  plane  seen  in  Figure  3-1.  The  X  axis 
is  chosen  to  coincide  with  the  ground  path  of  the  hypersonic  object  which  is  assumed  to 
moving  at  ground  speed,  v,  in  the  positive  X  direction.  Let  £z  and  £vbe  unit  vectors  in 

the  positive  X  and  Y  directions,  respectively.  Then,  V ,  the  horizontal  phase  velocity  of 
the  shock  front  is  given  by; 

y  =  v£z  1 

Now  let  t'  be  the  time  after  the  initial  contact  of  the  shock  front  with  the  surface  of  the 
earth  and  suppose  that  r0  identifies  a  point  which  moves  at  velocity,  V  ,  along  the  X  axis, 
i.e.; 

ra  =  vt'ez  2 

Assume  that  rs  locates  the  point  S(xs,  ys)  with  respect  to  a  moving  origin  at  r0 .  Then 

rs  =  +  y£y  3 


Now  let  r'  locate  the  same  point  with  respect  to  the  fixed  origin.  Then 


r'=r0+7s={xs  +  vt')ez  +  ysey  4 

Let  r  locate  the  observation  point,  0(x,  y)  with  respect  to  the  fixed  origin.  Then 

r=xsz  +  y£y  5 

Let  R  locate  the  O  with  respect  to  S.  Then, 

^  =  F_F'=(x-xf-vr')£J+(y-y,)ey  6 

Now  suppose  kc  is  a  horizontal  wave  number  vector  for  a  seismic  wave  which 
propagates  in  the  direction  <p  with  horizontal  phase  speed,  c.  Then, 

kc  =—  (cos<p£z  +sin<p£y}  7 

where  <p  is  measured  counterclockwise  with  respect  to  the  positive  X  axis;  (0  =  2nf  and  / 
is  the  frequency. 


The  Bow  Wave  Solution 

Let  P(ox  r  )  be  the  spectrum  of  the  overpressure  load  instantaneously  applied  at  the  point 
S  at  time  t'.  The  seismic  disturbance  created  by  this  load,  Su(R,t  -A,  as  seen  at  the  point 
O  at  some  later  time,  t  is  given  by 

Su(R,t- 1')  =  A  — ^ -  exp[-ia(t  -  «')]J exp(‘K  "  fyK  8 


where  A  is  a  constant  dependent  upon  the  elastic  properties  of  the  medium  and  q  is  the 
geometrical  spreading  coefficient.  Let  Au(R,t;T^jbe  the  sum  of  all  the  contributions 
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from  this  load  from  the  time  of  the  initial  contact  of  the  shock  front  with  the  surface  of 
the  earth  to  some  time,  T.  Then 

Au(r,,;T)  =  i;  aA^jLI  exp[-m(,  -  r')]J  exp(& „ .  9 

n  * 

Changing  the  order  of  integration  and  using  equations.  1-7  to  express  the  arguments  in 
equation  9  in  terms  of  their  component,  results  in 

Au(x  —  xs,y  —  ys,t;  vT)  =  AP(  co,  xs,ys)  exp\-i^-[vt  -{x-xs  )]Jtf(;c  -xs,y-ys;  vT) 
where 

(0 

^  r2ii  /— v- vj)  i*x-x. 

H(x-x„y-y,:vT)  =  l  e'  ■  '  - fta* 

[x!+(y-y,)2f 

It  is  shown  in  Part  II  that  for  large  vT  and  |r  —  <  vT ; 


J 

K 

(y-ys) 

cv 

H(x-xs,y-ys; vT.q)  =  B{x  -xs-vT,y-  ys;q)exp 
^{x ~xs’y-ys.t;vT)~  APB exp\-iax\exp\iKx (x - xs ) ± i Ky (y -  ys )] 


where 

CO 

K;=- 

C 


fc 

Vv 


,co 

K  =±  — 

>  c  1 


Hr 


10 


11 


12 


13 


14 


It  is  seen  from  equation  14  that  Au  can  take  one  of  two  different  forms  depending  on  the 
c  c 

value  of  — .  When  —  >7,  ktv is  imaginary.  Consequently,  as  illustrated  by  the  dashed 

line  in  Figure  3-2,  Au  becomes  an  elementary  inhomogeneous  wave  which  is  attached  to 
the  shock  front  at  point  (x„  yj.  It  propagates  with  speed,  v,  along  the  path  y  =ys  and  its 
amplitude  decays  exponentially  in  directions  normal  to  this  path.  On  the  other  hand,  if 
c 

—  <  1 ,  Kyis  real.  Therefore,  as  illustrated  by  the  solid  lines  in  Figure  3-2,  Au  describes 

two  elementary  plane  waves  which  are  attached  to  the  shock  front  at  point  (xs,  y5).  The 
elementary  plane  wavefronts  propagate  at  horizontal  speed,  c,  in  the  directions 
where 


15 
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Figure  3-2.  Diagram  of  the  elementary  seismic  bow  wavefronts  generated  at  the 
point  (Xsys  )  by  a  moving  shock  front.  Elementary  inhomogeneous  bow  waves 
(heavy  dashed  line)  are  generated  when  ov.  Elementary  propagating  plane  bow 
waves  (heavy  solid  lines)  are  generated  when  c<v. 

The  waves  generated  by  the  moving  point  pressure  load  are  referred  to  in  this  text  as 
elementary  seismic  "bow"  waves.  Now  let  u  be  the  resultant  of  the  elementary  bow 
waves  generated  at  all  the  of  the  points  on  the  shock  front  trace.  Now,  if  g(xs)  describes 
the  shock  front  trace  with  reference  to  an  origin  at  rs  and  ds  is  an  element  of  length  along 

the  shock  front  trace,  then; 

u(x,  y,Kx,Ky;  vT,  q)  =  AN(x,  y,  Kx,Ky ;  vT,  q,  g)  exp{-i(Ot)  exp[i(x:xx  +  Ky  y )]  16 

where 

N(x,y,Kz,Ky;vT,q)  =  J F(x-xs,y- g(x3);vT,q)exp[-i(Kxx1+Kyg(x3))\ds  17 
and 

F  =  P{(0,xs,g{xs))B(x-xs,y-g{xs);vT,q )  18 

Equation  16  is  the  formal  solution  for  the  seismic  response  to  a  moving  shock  front 
imposed  upon  the  surface  of  the  earth.  It  is  important  to  note  that  u  retains  the  phase 
behavior  of  Au,  that  is;  u  describes  two  propagating  plane  bow  waves  when  c<v  and  an 
inhomogeneous  bow  wave  when  c>v.  The  essential  difference  between  u  and  Au  is  the 
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occurrence  of  the 'complex  amplitude  factor,  N  on  the  right  hand  side  of  equation  16.  It 
can  be  seen  from  equation  18  that  F  is  a  complex,  non-linear  function  of  xs. 
Consequently,  it  would  be  difficult  to  obtain  closed  form  solution  for  N  given  reasonable 
assumptions  about  the  behavior  of  P  and  g.  Fortunately,  much  can  be  inferred  about  the 
behavior  of  N  by  the  examination  of  the  simple  model  of  the  shock  front  trace  which  is 
identified  below. 

As  shown  in  Figure  3-3,  let 

g(xs)  =  ±cot(d)(xs-d)  19 


Figure  3-3.  A  diagram  of  a  simplified  model  of  the  shock  front  trace. 


where  g(xs)  is  defined  only  in  the  interval,  b< xs<  d,  6  is  the  orientation  of  the  normal  to 
the  shock  front  trace,  and  the  point,  ( d,0 )  is  the  intersection  point  of  the  two  wings  of  the 
shock  front  trace.  Now  assume  that  c<  v  so  that  Ky  is  real.  Moreover,  for  the  sake  of 
simplicity  assume  that  P  is  independent  of  xs  and  that  as  shown  in  Figure  3-3,  the  shock 
front  trace  wings  are  of  length,  L.  Then,  neglecting  the  effects  of  geometric  spreading,  N 
for  this  particular  case  is  given  by; 

N  =  f  exP~  cot{e)(xs  -  <*))]*:,  20 

By  making  the  substitution  £=  xs-d  and  recognizing  that  b-d=-Lsin(d)  it  may  be  readily 
deduced  that; 
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21 


It  can  be  seen  by  inspection  of  equation  21  that 
\N\  _ 

pl  rr 


where 


22 


23 


which  illustrates  that  the  normalized  modulus  of  N  for  this  simple  model  of  the  shock 
front  trace  is  equal  to  the  modulus  of  a  sine  or  aperture  function.  The  normalized 
modulus  of  Nis  plotted  in  Figure  3-4  with  6  -  <pK  as  the  independent  variable. 


theta-phi  (degrees) 


Figure  3-4.  Normalized  Modulus  of  N  for  the  simplified  model  of  the  shock  front 
trace  (see  text  for  details). 

The  parameters  used  for  these  calculation  are  shown  in  the  panel  in  the  upper  left  comer 
of  the  figure.  They  are  appropriate  for  a  high  frequency  Pg  bow  wave  generated  by  a 
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shock  front  with  a  relatively  long  ground  path.  The  results  seen  in  Figure  3-4  demonstrate 
that  N  is  maximized  when  the  normal  to  the  trace  of  the  shock  wave  front  is  closely 
aligned  with  k  and  diminishes  rapidly  as  the  misalignment  increases.  In  particular,  it 
should  be  noted  that  for  this  set  of  model  parameters,  that  if  the  misalignment  is  greater 
than  about  ±4° ,  then  the  normalized  value  of  |A/j  is  down  more  than  25  dB  with  respect 

to  its  peak  value.  These  results  indicate  that  the  elementary  bow  waves  generated  by  the 
point  sources  distributed  along  the  shock  front  trace  will  be  strongly  attenuated  by 
destructive  interference  unless  the  normal  to  the  shock  front  traces  happens  to  be 
approximately  parallel  to  k  . 

A  similar  behavior  can  be  reasonably  anticipated  that  for  actual  shock  front  traces. 
However,  in  .this  case  Q  is  not  a  constant,  but  rather  it  is  a  continuous  variable  defined  in 
the  range: 

-  cos'1  \  —  |  <6<  cos'1  f  —  ]  23 

V  v  J  {  v  J 

where  ca  is  the  local  sound  speed.  Consequently,  it  is  argued  that  for  actual  shock  fronts 
with  relatively  long  ground  paths,  destructive  interference  will  strongly  attenuate  the 
elementary  bow  waves  generated  everywhere  along  the  trace  except  at  those  points 
where  6  =  (pK.  The  spatial  relationship  between  the  shock  front  and  its  associated  bow 
wave  implied  by  this  argument  are  conceptually  illustrated  in  Figure  3-5. 


Figure  3-5.  Spatial  relationship  between  the  shock  front  trace  (dashed)  and  its 
associated  bow  wave  (solid). 


2 


It  can  be  seen  from  this  figure  that  the  bow  wave  front  is  tangent  to  the  shock  front  at  the 
points  where  the  6  =  —-0K  and  intersects  the  X  axis  at  the  point  (d,0).  It  is  also  important 

to  note  that  if  c(6)  denotes  the  speed  of  the  shock  front  in  the  direction  6,  then 
c{6 )  =  v  cos(d)  24 

and  at  the  points  where  6  =  ±<pK 

c(6)  =  v  cos(±<p  K)  =  c  25 

In  other  words,  when  c<  v,  bow  waves  will  "couple"  to  the  shock  front  at  those  points 
where  its  horizontal  phase  speed  is  equal  to  the  horizontal  speed  of  seismic  waves  in  the 
underlying  medium.  Finally,  based  upon  this  argument,  it  follows  that 

u  ~  AW exp\iKx{x  -  J)]exp[iJCyy]  26 

where  W  is  a  function  which  accounts  for  the  overpressure  distribution  along  the  shock 
front  trace  as  well  as  the  effects  of  geometrical  spreading. 
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PART  n.  SEISMIC  "BOW"  WAVES  GENERATED  BY  MOVING  POINT 
PRESSURE  LOADS  APPLIED  AT  THE  SURFACE  OF  THE  EARTH 


It  is  shown  in  the  accompanying  paper  that  the  earth  motion  observed  at  time,  t,  and  at  the 
point  {x,y)  due  to  a  point  pressure  load  which  has  moved  at  speed  v  on  the  line  y=ys 
parallel  to  the  X  axis  over  a  period  of  time  T  may  be  expressed  as 


Au  =  A P(a>)exp<  -ico  t 


1 


where 


fieos*^)x 

■ - jd%dtp  2 

x2+{y-ys)2]2 

In  order  to  simplify  the  notation,  let 
xs=ys=0 


H= 

JO  Jx-x. 


-vT  , 


kx  =  — cosip 
c 

7  ® 

kr  =  — simp 
c 


L  —  vT 


Now  let  Wbe  defined  by; 


3 


then 

JO 


exp\ 


f 

\ 

—  - 

V  2 

YVC  ) 

m 


6) 


-kx;x,y,L  \dkz 


4 


It  can  be  seen  from  equation  4  that  H  may  be  identified  as  a  convolution  in  the  wave 
number  domain  between  a  plane  wave  and  the  function,  W.  Moreover,  it  may  be  seen 
from  equation  3  that  W  may  be  viewed  as  the  wave  number  spectrum  of  the  output  of  a 
phased  linear  source  array  of  length,  L,  whose  individual  inputs  are  amplitude  weighted 
by  the  geometrical  spreading  factor.  An  expression  for  W  is  derived  in  the  following 
section. 

Evaluation  of  W  for  Body  Wave  Inputs 

Let 


Then,  it  follows  directly  from  equation  3  that 

where 

f\ = Lexp(-^yz 


5 

6 

7 
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is  the  wave  number  spectrum  of  a  linear  source  array,  and 
.  _  r  exp  (-ha). 


(x!+?T 

is  the  one  dimensional  wave  number  spectrum  of  the  geometrical  spreading  factor.  It  is 
straightforward  to  show  that: 


•  frjL 

,  „  .  (  L\\Sin{  2 

/1=2exp^-IJj— — 

and  assuming 


which  is  valid  for  far  field  body  waves  generated  by  a  surface  point  source,  the  calculus 
of  residues  may  be  used  to  show  that: 

exp(-|%)  w 


IV  =  ij%xp  sm^-yjexpt-liJ-Aly)- 


(  XL i 

By  replacing  sin  —  in  equation  1 1  with  its  exponential  representation,  it  can  be  shown 

2  J 


W  =  -^{exp(77y)[£1(z1*T7)  -  Ei(zlv)]  ~  exp^T^E,  (-z,r?)  -  Ex  (-z2r?)]} 


where 


is  the  exponential  integral;  and 

Zl=3,-“  ,  14 

z2  =  y-i(x-L) 

and  the*  indicates  the  complex  conjugate. 

For  small  values  of  r\,Ej  may  be  safely  represented  by  a  Maclaurin  series  of  the  form 
£,(f)  =  -D- log,(?)-jT-£-  15 

n=l  n'n‘ 

where  D  is  the  Euler-Mascheroni  constant.  Substitution  of  this  series  into  equation  12 
yields; 


Re  W  =  +  smh(w)±  -  cosh^gig^L- 


!m  W  =  i  sinhMa  -  cosh(rjy)±  ( j " 
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where 


cc  =  Re  log. 


P  =  ]m  log  -2- 


r«,=Re[z,"]-Re[Zl-] 
imz,  =Im[z,,']-Im[z!'] 


It  can  be  seen  from  these  equations  that 
w(-n)-w(-n) 

For  larger  values  of  jj,  Ei  should  be  represented  by  its  asymptotic  expansion 

£,(f) .  exPK)[i-^+|y-|l+-(-r  <^| 


In  this  case  substitution  of  the  asymptotic  expansion  into  equation  12  yields 

y{  n  t\2  t?  K  >  vn  l 

hnW  =  +  2  iilnl  lV'1 

y{  V  V2  rf  v  '  77"  J 

where 

cn  (77)  =  sin(^x)  Re(z,"" )  -  sin[77(x  -  L)]  Re(z ;n ) ;  n  odd 
cn  (77)  =  cos(77x)lm(z“" )  -  cos[77(x  -  l)]lm(zj" ) ;  n  even 
Znil)  =  cos(77x)Re(z1"")  -  cos[77(x  -  L)]Re(z2"n) ;  n  odd 

=  sin(77y)lm(z-n)  -sin[77(x-  L)]lm(z2_n) ;  n  even 

The  behavior  of  W  for  large  L 

For  L>r  and  for  77  near  zero,  it  is  found  from  equations  16  that 

ReW- - L--Lrf  -0(l3t]4)— ■ 

r  sind  2  '  ' 

ImW  =  Cf7j-o(L2773) — 

where 

e-um"U' 1 


r  =  (x2  +  y2)? 

For  larger  values  of  77,  it  is  found  from  equations  20  and  21  that 


76 


’  sin  (777- cos  0)  sin  0  cos(77rcos0)sin20  „  sin(77rcos0)sin  30 


Re  W 


rsin0  f  r2 

...A- 


’  cos(7jr  cos  0)  sin  0  sin(77r  cos  0)  sin  20  cosjrjrcos  0)  sin  30 


ImW  = 


rsin0 


'§ 


It  can  be  seen  from  these  equations  that  |W(t7)|  will  have  an  absolute  maximum  at  77  =  0 

and  will  diminish  as  -Lrf  for  77  near  zero,  and  at  least  as  rapidly  as  (rT])'1  for  larger 
values  of  77.  Thus  for  large  L  and  L>r,  W  may  be  identified  as  a  beam  function  whose 
half  width,  rjw  is  dependent  on  L  and  r . 

According  to  convention,  rjw  is  the  smallest  real  root  of  : 

|wf  =inm(|wf)  25 

Let  <50*  be  the  beam  width  expressed  in  degrees.  Then 

W<>)-0(n.))  26 


cos[«(j))]  =  c^--| 
Thus  at  77  =  0 


cos[0(O)]  =  — 
while  at  77  =  77  K 


Therefore,  360  _,([  1  Vw 

56  = -  cos - cos  c - 

K  v  v  J  Uv  (0 


Equations  28,  29  and  30  may  be  used  to  derive  approximate  values  for  <50*.  An  example 
of  these  derivations  is  shown  in  Figure  3-6.  The  curves  shown  in  this  figure  illustrate  the 
influence  of  increasing  values  of  r  and  L  on  5<pw  for  a  model  which  is  appropriate  for  a 
high  speed  source  path  near  a  station  located  at  the  surface  of  a  medium  with  a  P  or  S 
wave  velocity  of  4.5  km/sec. 
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Figure  3-6.  The  beam  width  of  W  as  function  of  r  and  L. 


i  i  ^ 

Notice  that  8<pw  decreases  as  —  for  L<r  and  as  ~T£\  f°r  L>  r-  It  is  especially 

significant  to  note  that  for  r>1000  km.,  5^  is  less  than  about  0.2  degree.  Clearly  then, 
for  large  r  and  L>r,W  will  define  a  sharply  tuned  beam  focused  at  77  =  0.  This  behavior 
is  illustrated  in  greater  detail  by  the  results  shown  in  Figure  3-7.  The  upper  panel  of  this 
figure  portrays  the  modulus  of  W(rj)  calculated  for  a  model  with  the  parameters  listed  in 
the  left  hand  side  of  the  panel.  The  phase  of  W(ti)  is  plotted  in  the  lower  panel.  In  order 
to  rapidly  evaluate  W  for  a  wide  range  of  77,  the  Fast  Fourier  Transform  (FFT)  algorithm 
was  applied  equation  12  to  produce  the  data  shown  in  this  figure.  To  facilitate 
comparisons  between  different  wave  types,  the  modulus  has  normalized  by  dividing  by 
through  by  its  maximum  value.  It  is  seen  from  equation  22-a  that  for  a  body  wave  this 
.  .  B 

factor  is  — — - .  As  noted  earlier,  a  spreading  coefficient  of  2  is  appropriate  for  a  far 
r  sin#  r 

field  body  wave  from  a  surface  point  source. 

Observe  that  for  this  case 

\W (7j)|  <  0.1  max(|W (77)]) ;  |t7|  >  0.005  31 

which  indicates  that  it  is  a  very  sharply  tuned  beam  function. 
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Example  of  W  for  a  Surface  Wave  Input 

The  beam  function  for  a  far  field  surface  wave  from  a  surface  point  source  is  shown  in 
Figure  3-8.  It  was  generated  by  changing  the  spreading  coefficient,  q,  from  2  to  0.5  and 
reapplying  the  FFT  algorithm  to  equation  3,  using  the  same  model  parameters. 

Notice  that  the  width  of  the  modulus  of  the  surface  wave  beam  function  is  significandy 
smaller  than  the  body  wave  beam  width.  This  reflects  the  fact  that  surface  wave 
geometrical  spreading  is  substantially  less  severe  than  body  wave  geometrical  spreading. 


Figure  3-7.  An  example  of  the  Beam  Function  for  body  waves  generated  by  a 
moving  surface  point  pressure  load. 


Figure  3-8  .  An  example  of  the  Beam  Function  for  surface  waves  generated  by  a 
moving  surface  point  pressure  load. 
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An  Approximate  Expression  for  H  for  Body  Wave  Inputs 

Based  upon  the  evidence  and  arguments  presented  previously,  it  may  be  concluded  that 
for  large  r  and  L>  r,  both  body  and  surface  wave  beam  functions  may  be  approximated  by 
functions  of  the  form; 


W(ri)  =  B(x,y;L,q )  77  =  0 


where  in  each  case  B  is  the  value  of  |W|  at  77=0.  From  equations  23,  it  can  be  seen  that 
for  body  waves,  B  is  given  by 


Consequently,  substitution  of  equations  32  and  33  into  equation  4  and  performing  the 
indicated  convolution  results  in  the  following  expression  for  U  created  by  body  waves 


34 


Equation  34  exhibits  the  essential  properties  of  the  seismic  "bow  wave”  generated  by  a 
moving  pressure  load.  Observe  that  if  v  >  c  the  radical  in  equation  42  is  real  and  H 
becomes  a  plane  wave  which  propagates  at  speed  c  in  the  direction 


<f>  =  cos 
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On  the  other  hand,  if  v  <  c ,  the  radical  in  equation  34  is  a  pure  imaginary  and  H 
becomes  a  disturbance  which  propagates  at  the  speed  of,  and  on  the  path  of  the  moving 
load  but  decays  exponentially  in  a  direction  normal  to  the  path. 

The  evolution  of  H  as  L  increases  is  illustrated  by  the  relative  amplitude  profiles  shown 
in  Figure  3-9-a.  The  relative  amplitude.  A,  of  His  defined  as 
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which  is  seen  to  be  the  ratio  of  the  amplitude  of  H  to  the  amplitude  of  a  far  field  body 
wave  generated  by  the  point  pressure  load  at  the  origin.  The  profiles  were  calculated 
along  a  line  parallel  to  the  load  path  but  offset  100  km  in  order  to  approximate  far  field 
conditions.  Notice  that  for  distances  along  the  profile  which  are  much  greater  than  the 
path  length,  the  amplitude  of  H  scales  as  the  amplitude  of  a  point  source  body  wave  of 
strength  L  ;  whereas  for  distances  comparable  to  or  less  than  the  path  length  the 
amplitude  of  H  differs  significantly  from  that  of  a  point  source  body  wave  particularly 
for  large  L.  The  two  dimensional  behavior  of  the  relative  amplitude  of  H  for  large  L  is 
was  also  studied.  The  data  in  this  case  were  calculated  for  the  case  where  L=  1000  km. 
Inside  the  region,  r<L,  the  relative  amplitude  of  H  reaches  a  maximum  of  about  85  dB  for 
x  near  L  and  y=100  km.,  then  approaches  0  dB  as  x  and  y  approach  zero. 
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Figure  3-9.  Relative  amplitudes  of  H.  Relative  amplitude  profiles  of  U  taken 
parallel  to  the  load  path  at  an  offset  distance  of  100  km. 


Finally,  by  letting 
v  =  x-L 

and  recognizing  that 
L  —  vt 
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equation  34  may  rewritten  as 
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which  for  large  values  of  t  reduces  to 
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Equation  40  demonstrates  that  for  large  values  of  L  and  as  seen  in  a  reference  system 
fixed  to  the  moving  load  ,  H  is  essentially  independent  of  time. 

It  may  be  concluded  from  this  analysis  that  relatively  weak  pressure  sources  will  create 
signals  that  are  detectable  at  stations  which  are  located  at  within  a  few  hundred 
kilometers  of  the  load  path  when  the  length  of  the  path  is  of  the  order  of  1000  km  and 


v  >  c . 
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APPENDIX4.  Preliminary  Tests  of  Alternative  Infrasonic 

Sensor  Designs 

Introduction 

As  a  result  of  the  Comprehensive  Test  Ban  Treaty  (CTBT),  there  is  a  renewed 
interest  in  detecting  and  recording  infrasonic  phenomena.  Up  to  now,  only  a  few  research 
and  volcano  observatories  have  maintained  infrasound  sensors  and  data  collection  facilities. 
These  existing  stations  have  not  generally  recorded  high  dynamic  range  (118  dB), 
broadband  (0.02-16  Hz)  signals.  Experimental  infrasound  arrays  have  been  deployed,  but 
these  do  not  always  have  the  CTBT  recommended  1-2  km  aperture. 

Current  operational  systems  use  one  of  two  types  of  microbarographs,  the  absolute 
pressure  device  (Figure  1)  and  the  differential  pressure  microbarograph  (Figure  2)  used  at 
Lajitas,  Texas  (TXAR),  and  Warramunga,  Australia  (WRA).  Each  has  advantages.  The 
differential  microbarograph  is  easy  to  assemble  from  off-the-shelf  commercial  products, 
making  it  inexpensive;  measures  only  the  band  of  interest,  making  it  suitable  for  low 
dynamic  range  recorders;  is  capable  of  a  bandwidths  that  include  audio,  allowing  overlap 
calibration  with  microphones;  and  is  reasonably  linear  over  wide  variations  in  atmospheric 
pressure,  minimizing  the  need  to  record  atmospheric  pressure.  The  absolute 
microbarograph  is  much  less  sensitive  to  ambient  temperature  and  temperature  gradients, 
allowing  it  to  be  used  in  shallow  vaults;  delivers  state-of-health  atmospheric  barometric 
pressure  as  part  of  the  data;  and  is  simpler  to  describe,  giving  a  simplified  theoretical 
response. 

Work  at  SMU  has  concentrated  exclusively  on  differential  sensors  because  they  are 
reasonably  priced,  are  physically  compact,  and  are  compatiable  with  existing  seismic 
digitizers.  The  microbarographs  operating  at  TXAR  were  assembled  for  nominal  material 
cost  using  simple  tools.  The  compact  assembled  system  is  reasonably  rugged  and  fits 
inside  the  8  in.  borehole  above  the  seismometer.  Since  the  electronics  and  interface  are 
plug  compatible  with  existing  seismic  digitizers,  a  seismic  array  may  be  easily  upgraded  to 
a  seismo-acoustic  array  by  replacing  a  single  component  digitizer  with  a  two  channel 
digitzer. 

Previous  Work 

Specific  work  in  practical  infrasonic  instrumentation  design  is  sparse.  Original 
theoretical  work  by  Benade  (1968),  Burridge  (1971),  Cook  and  Bedard  (1971),  Daniels 
(1950  and  1959)  describe  the  basic  physics  involved  in  microbarograph  responses  and  the 
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responses  of  the  associated  noise  reducing  pipe  arrays.  The  work  by  Haak  and  Wilde 
(1996)  is  an  unusually  complete  summary  of  the  physics  for  microbarograph  design. 

Like  much  of  the  infrasound  community,  SMU  has  used  porous  garden  soaker 
hoses  in  noise  reducing  pipe  arrays  as  an  inexpensive  alternative  to  the  more  elaborate 
Daniel’s  line  microphone.  A  previous  study  at  SMU,  (Hayward,  1996)  considered  the 
characteristics  of  the  sensors  and  hoses  currently  deployed  at  TXAR  and  discussed  the 
response  variability  and  resonance  peaks  as  a  function  of  differing  hoses  (all  of  the  same 
manufacture).  The  study  also  considered  several  methods  for  characterizing  the  response 
of  a  porous  hose.  For  example,  responses  of  two  identical  hoses  laid  outdoors  parallel  to 
each  other  7  cm  apart  was  incoherent  over  much  of  the  spectrum.  A  differential  sensor  was 
used  to  measure  the  pressure  gradient  along  the  hose  and  this  was  compared  against  the 
readings  from  a  differential  microbarograph.  The  two  relation  between  the  two  signals  was 
complex  and  unexplained.  Finally,  sensor  noise  levels  and  responses  measured  for  three 
microbarographs  subsequently  deployed  at  TXAR  suggested  that  noise  levels  for  the 
sensors  could  vary  widely.  Resolving  these  problems  was  postponed  until  the 
characteristics  of  the  equipment  and  processes  were  understood. 

Since  the  initial  experiments  with  sensors  at  SMU,  Haak  and  Wilde  (1996) 
reviewed  the  physics  involved  in  porous  hose  design  and  suggested  alternative  designs 
such  as  those  based  on  a  closed  compliant  tube,  and  McLaughlin  (1997)  analyzed  various 
porous  hose  effects  based  on  simple  laboratory  measurements  and  the  results  from  a  finite 
difference  modeling  program. 

Current  Sensor  Design 

The  microbarograph  currently  used  at  TXAR  (Figure  3)  is  based  on  Reinke’s 
(1985)  design  which  in  turn  was  based  on  the  early  work  of  Stachura,  Siskind,  and 
Engeler  (198 1)  utilizing  differential  pressure  gauges  to  measure  blasting  overpressures.  A 
short  section  of  iron  pipe  forms  the  backing  volume.  Instead  of  a  capillary  tube  as  a  bleed 
resistence,  a  variable  needle  valve  is  used  to  set  the  low  frequency  comer  of  1/60  Hz.  The 
valve  is  adjusted  in  the  lab  and  then  locked  into  place  with  a  small  dab  of  glue.  A  standard 
garden  hose  connector  at  the  top  allows  connection  to  the  pipe  array.  At  TXAR,  the 
sensors  are  placed  about  2  m  below  the  surface  for  security  and  temperature  stability. 

The  pressure  sensor  is  the  Validyne  P305D  sensor,  a  differential  pressure  gauge 
with  a  range  of  ±  0.125  PSI  (±850  Pa)  and  a  ±5  Volt  output.  Higher  pressure  ranges  are 
available,  but  this  would  limit  the  use  for  small  signals.  Careful  handling  during  installation 
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is  required  to  avoid  exceeding  the  rated  pressure  of  ±3  inches  water.  While  overpressures 
generally  do  no  permanent  harm,  they  can  disturb  the  calibration. 

The  sensor  itself  is  a  variable  reluctance  device.  Two  coils  are  positioned  on  either 
side  of  a  magnetically  permeable  stainless  steel  diaphragm.  Small  changes  in  the 
diaphragm  position  cause  a  change  in  the  inductance  of  the  sensor  coils,  an  increase  on  one 
side  and  a  decrease  on  the  opposite  side.  There  are  a  number  of  methods  to  convert  this 
impedance  change  to  an  analog  signal.  Most  techniques  incorporate  the  two  sensor  coils 
into  an  AC  bridge  balanced  at  zero  pressure.  The  P305  electronics  uses  an  internal  5.2 

KHz  ±5  Volt  square  wave  excitation  and  produces  an  output  with  less  than  -70  dBV  (0.05 
Pa  equivalent)  noise  at  the  half  the  primaiy  excitation  frequency  (2600  Hz).  Antialias 
filtering  in  the  digitizer  removes  the  remnants.  The  gauges  produces  a  ±5  Volt  single  ended 

signal  referenced  to  the  negative  supply  bandlimited  from  DC  to  200  Hz  from  an  unfiltered 
primary  12  Volt  power  supply.  Sensor  output  is  digitized  at  40  SPS  on  a  24-bit  Science 
Horizons  digitizer  with  a  2  fiV  (0.00035  Pa)  LSB. 

The  construction  and  plumbing  in  the  current  TXAR  design  was  based  on 
conservative  modification  and  quick  fabrication  of  Reinke’s  design  rather  than  by  extensive 

calculation  and  testing.  It  resulted  in  a  system  that  possessed  a  number  of  unknowns  and 
idiosyncrasies.  In  spite  of  this,  for  the  past  two  years  TXAR  has  produced  highly  relevant 
infrasonic  data  leading  to  significant  advances  in  the  understanding  of  infrasonic  processes 
above  1/2  Hz. 

Problems  with  Current  System  at  TXAR 

The  microbarograph  design  can  be  divided  into  three  logical  subsystems,  the  sensor 
and  associated  signal  processing  electronics,  the  backing  volume  and  associated  plumbing 
required  to  create  a  low  cut  filter  (the  microbarograph  itself),  and  the  noise  reducing  pipe 
array.  While  all  of  these  components  contribute  to  the  total  system  response,  most  obvious 
problems  in  the  current  design  can  be  attributed  to  one  subcomponent. 

Sensor 

There  are  four  problems  with  the  Validyne  P305  sensor;  high  self  noise,  DC  drifts, 
single  ended  outputs,  and  a  strong  sensitivity  to  RF  interference.  None  of  these  preclude 
the  use  of  the  sensor,  but  all  have  to  be  considered  and  handled  during  testing  and 
installation. 
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During  the  initial  noise  tests,  one  sensor  was  found  to  have  a  noise  floor  14  dB 
higher  than  the  quietest.  At  the  time,  this  was  attributed  to  a  noisy  set  of  electronics,  and  it 
was  felt  that  it  might  be  necessary  to  select  units  prior  to  deployment.  The  noisy  unit  was 
subsequently  retained  as  a  laboratory  test  unit.  The  manufacture’s  specifications  suggested 

only  a  RMS  noise  floor  <5  mV  (equivalent  to  0.86  Pa).  Measured  noise  in  the  band  of 
interest  (0.02-20  Hz)  was  far  below  this. 

Sensor  output  is  relative  to  the  negative  power  supply.  In  the  current  AIM  systems, 
the  primary  power  supply  is  completely  isolated  from  signal  paths  by  DC-DC  converters. 
When  a  Validyne  sensor  is  installed  in  the  system  though,  the  negative  primary  supply 
must  be  tied  to  negative  signal  output.  This  results  in  a  signal  path  back  through  the  power 
supply  that  can  introduce  noise  into  the  system. 

When  small  amounts  of  modulated  RF  energy  above  100  MHz  are  present  in  the 
cable  connecting  the  Validyne  electronics,  the  output  may  be  excessively  noisy.  This  was 
noted  at  TXAR  for  sites  using  packetized  RF  telemetry.  The  effect  does  not  seem  to  be 
consistent.  The  digitizer  front  ends  also  are  susceptibile  making  it  difficult  to  isolate  the 
problem.  In  the  field,  the  RF  interference  is  minimized  by  moving  antennas  away  from  the 
electronics,  shielding  the  borehole,  and  nulling  the  local  field  by  adjusting  cable  positions. 

DC  offset  varies  with  ambient  temperature  partly  from  offsets  in  the  electronics  and 
partly  due  to  minor  changes  in  the  diaphragm  position  as  the  sensor  block  is  subjected  to 
differing  temperature  gradients.  Offset  also  varies  with  sensor  orientation,  since  gravity 
affects  diaphragm  position.  This  DC  offset  may  be  set  to  zero  with  an  internal  adjustment; 
however,  in  outdoor  field  experiments  with  the  sensor  on  the  surface,  it  was  not  possible  to 
keep  the  offset  within  62  mV  over  a  several  hour  period.  This  offset  is  a  problem  for  16- 
bit  digitizers,  which  have  a  maximum  input  of  62  mV  (32767  *  2  uV)  when  the  sensitivity 
is  low  enough  to  detect  small  signals.  For  laboratory  and  field  testing,  sensors  were 
adjusted  for  zero  offset  and  tests  were  kept  brief.  When  installed  at  TXAR,  the  systems 
were  connected  to  24-bit  digitizers. 

Microbarograph 

Problems  with  the  microbarograph  are  divided  into  three  groups;  unknown 
characteristics,  practical  problems,  and  technical  problems  with  the  response. 
Characteristics  such  as  impedance,  thermal  constants,  backing  chamber  volume  and  fore 
volume,  are  unknown.  Practical  problems  include  the  fragile  nature  of  the  finished 
product,  the  tendency  of  the  adjustable  needle  valve  to  get  bumped,  and  a  difficulty  in 
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keeping  critical  areas  of  the  plumbing  clear  of  dirt  and  moisture.  There  are  three  specific 
technical  problems  with  the  response;  the  construction  results  in  a  high-cut  acoustic  filter  at 
about  5  Hz  (although  it  varies  slightly  from  unit  to  unit),  units  in  the  field  drift  from  the 
laboratory  setting  and  eventually  have  wildly  different  responses,  and  there  is  no  way  to 
remotely  calibrate  the  unit. 

It  is  essential  to  keep  the  plumbing  clean  and  dry.  Once  assembled,  but  prior  to 
connecting  the  sensor,  the  plumbing  is  washed  with  a  degreasing  cleaner  to  eliminate  all 
traces  of  grease  and  dirt.  Joints  are  sealed  with  Teflon  tape.  On  repeated  assembly  and 
disassembly  small  pieces  of  tape  have  a  tendency  to  collect  inside  fittings. 

The  needle  valve,  used  to  adjust  the  response,  is  particularly  troublesome.  As 
delivered,  the  spring  loaded  needle  is  greased  at  several  points.  In  order  to  obtain  the 
proper  response  the  valve  had  to  be  nearly  closed.  The  resulting  orifice  is  inconsistent. 
Grease  in  the  spring  assembly  caused  valve  creep  over  periods  of  minutes  and  hours.  It  is 
essential  to  completely  disasemble  the  needle  valve  and  throughly  degrease  the 
components.  Still,  valves  set  correctly  in  the  morning  can  creep  completely  closed  by  the 
evening.  Mechanical  backlash  in  the  system  requires  that  the  adjustments  be  made  from 
one  direction  (open  to  closed). 

Once  the  valve  is  properly  set,  the  control  is  cemented  in  place.  This  did  not  always 
result  in  a  fixed  setting.  Moderate  blows  to  the  assembly  or  large  temperature  changes 
sometimes  resulted  in  a  significant  change  in  the  response 

Laboratory  tests  demonstrated  that  the  system  has  an  unplanned  5  Hz  low  pass 
acoustic  filter.  This  was  demonstrated  by  comparing  the  response  of  a  plumbed  system 
with  a  sensor  running  as  a  simple  microphone  (using  a  closed  backing  volume)  to  noise 
generated  by  a  high  pressure  air  jet.  Although  few  infrasonic  signals  of  interest  have 
significant  energy  above  5  Hz  and  the  addition  of  noise  reducing  pipes  will  further  attenuate 
these  higher  frequencies,  this  additional  filter  term  complicates  the  system  response. 

To  a  first  order  the  response  of  the  microbarograph  is  a  simple  1  pole  RC  filter.  It 
was  intended  that  the  comer  be  at  1/60  Hz.  Observation  of  field  data  suggests  that  two  of 
the  sensors  now  have  comers  at  far  longer  periods.  Initially  this  appeared  as  a  slow  drift 
towards  longer  frequencies.  It  is  unknown  if  the  present  position  is  stable. 

Sensors  have  never  been  calibrated  for  sensitivity  other  than  the  initial  factory 
calibration.  Factory  calibration  of  ±0.25%  is  sufficient  given  the  uncertainties  in  the  system 
response,  provided  no  damage  or  overpressure  has  occurred.  Overpressures  of  as  little  as 
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10  times  (18  KPa)  may  change  the  sensitivity  by  2%.  Pressures  of  100  times  (180  KPa) 
may  make  a  4%  change.  A  more  critical  problem  is  dirt  or  moisture  in  the  diaphragm 
chamber  which  may  result  in  sensitivity  changes  of  100%  or  more.  Being  able  to  remotely 
calibrate  sensors  would  significantly  increase  the  confidence  in  the  response  and 
sensitivity. 

The  acoustic  and  physical  characteristics  of  the  plumbing  are  unknown  and 
dimensions  that  could  be  used  to  estimate  these  are  unknown  and  uncontrolled.  While 
simple  experiments  were  done  to  measure  the  sensor  noise  characteristics  and  dynamic 
response  of  the  system  to  large  signals,  more  specific  measurements,  such  as  acoustic 
impedances,  thermal  capacity,  or  internal  resonanaces  were  undetermined.  Without  these, 
designing  a  matching  feed  line  for  a  pipe  array  is  difficult. 

Pipe  Array 

The  pipe  array  at  TXAR  consists  of  five  or  six1  5/8  in.  by  50  ft.,  porous  garden 
hoses  connected  to  25  ft.  of  5/8  in.  hard  wall  garden  hose  in  turn  connected  to  a  manifold 
with  a  volume  of  roughly  2  liters.  A  25  ft.  hard  wall  garden  hose  connects  the  manifold  to 
the  microbarograph  hanging  2m.  below  the  surface  in  an  8  in.  cased  borehole  shared  with  a 
seismometer.  Like  the  microbarograph,  the  pipe  array  was  constructed  with  easily  obtained 
local  material  and  deployed  quickly.  Based  on  initial  tests  on  the  SMU  campus  and 
observations  of  data  collected  at  TXAR,  the  array  has  the  following  problems:  strong 
variable  hose  resonances,  incoherence  on  adjacent  hoses,  and  leaks  of  end  points.  In 
addition  the  hose  state-of-health,  acoustic  characteristics,  spatial  and  frequency  response, 
and  aging  characteristics  are  unknown. 

At  TXAR  sites,  pipe  arrays  appear  to  have  unique  resonances,  although  all  sites  use 
the  same  array  design.  The  origin  resonances  have  not  been  investigated. 

Hose  connection  points  often  have  substantail  air  leaks  as  the  result  of  a  poorly 
seated  washer  or  an  incomplete  crimp  on  the  connector  ferrule.  Some  leaks  were 
equivalent  to  the  contribution  of  50  ft.  of  hose.  Leaks  were  located  by  observing  the  signal 
from  a  low  pressure  air  jet  as  it  was  moved  along  the  hose  and  connector.  Loose 
connections  were  reseated  and  tightened  while  leaky  ferrules  were  taped. 

Soaker  hoses  at  TXAR  are  subject  to  intense  sunlight,  dust,  and  occasional 
investigation  by  javalena  (wild  pigs)  in  the  area.  In  addition  aggressive  vegetation  and 
insects  are  risks.  Because  the  acoustic  resistance  of  the  pipe  array  is  significant,  a  single 

1  One  site  uses  six  hoses,  the  other  three  use  five  hoses. 
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cut  or  tear  in  any  hose  will  result  in  a  change  in  the  frequency  response  and  a  significant 
reduction  in  the  spatial  filtering. 

Unfortunately,  it  is  nearly  impossible  to  monitor  the  state  of  the  hose  array  from 
observing  the  data  itself.  If  there  were  archived  background  responses  during  no  wind 
conditions,  it  might  be  possible  to  recognize  a  change  in  the  spectra  that  would  indicate  a 
cut  or  break  in  the  hose,  but  wind  conditions  are  not  recorded  nor  has  noise  spectra  been 
routinely  collected.  Examining  the  hoses  on  a  routine  maintenance  visit  could  locate 
obvious  tears  and  breaks,  but  could  also  miss  small  punctures  and  would  require 
significant  time.  It  would  probably  be  more  economical  to  simply  replace  hoses — if  one 
could  be  assured  the  acoustic  response  of  the  system  would  be  the  same. 

Experiments  by  Haak  and  Wilde  (1996)  and  McLaughlin  (1997)  have  confirmed 
that  acoustic  characteristics  of  soaker  hoses  change  substantially  as  a  result  of  aging  and 
moisture  in  exposed  hoses.  This  would  affect  both  the  frequency  and  spatial  response  of  a 
system.  Haak  and  Wilde  (1996)  suggested  that  hoses  be  protected  from  moisture  and  dirt 
by  elevating  them  from  the  ground  and  adding  a  roof.  This  though  has  the  disadvantage  of 
moving  them  farther  from  the  ground  boundary  layer. 

Purpose 

The  purpose  of  this  paper  is  to  document  an  improved  microbarograph  design,  to 
discuss  the  initial  laboratory  experiments  using  the  new  design  and  to  suggest  methods  for 
field  calibration  and  state-of-health  monitoring. 

The  improved  design  involves  different  plumbing  of  the  same  P305D  sensor,  and 
an  initial  consideration  of  two  alternative  sensors,  one  with  an  enhanced  sensitivity  (up  to 
30  times)  and  one  with  a  lower  price  (about  1/5). 

This  is  intended  as  a  starting  point  for  further  modifications,  for  comparing 
alternative  sensors,  for  measuring  the  acoustic  and  electrical  characteristics  of  the  system, 
and  for  comparing  the  measured  responses  with  theoretical  responses. 

Three  simple  calibrators  for  remote  field  calibration  and  state-of-health  monitoring 
are  suggested  and  demonstrated  on  the  prototype  system. 
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Results 


Microbarograph  Design 

The  design  goals  were  to  increase  the  backing  volume  to  1  liter;  to  replace  the 
variable  needle  valve  with  either  a  fixed  orifice  or  capillary  tube;  to  allow  some  form  of 
field  calibration;  to  increase  the  high  frequency  cutoff  of  the  system  beyond  the  digitizer 
cutoff  of  20  Hz;  to  increase  the  thermal  mass  such  that  the  sensor  is  less  sensitive  to  small 
temperature  changes;  to  keep  the  system  small  enough  to  fit  within  the  existing  boreholes  at 
TXAR  (8”  internal  diameter  with  sufficient  room  to  clear  existing  cables);  and  to  keep  the 
design  flexible  and  inexpensive  such  that  changes  could  be  made  during  tests.  Like  the 
TXAR  design,  this  improvement  uses  the  Validyne  P305D  sensor. 

The  design  is  based  on  4  in.  schedule  40  PVC  tubing  (Figures  4,  5,  6,  and  7). 

PVC  prototypes  may  be  quickly  constructed,  are  inexpensive,  and  provide  a  significant 
thermal  insulation  between  the  outside  and  inside  air  temperatures.  The  inside  of  4  in.  PVC 
end  caps  have  been  reamed  slightly  to  close  each  end  of  a  tube  with  a  slip  fit.  O-rings  in 
machined  grooves  on  the  body  of  the  instrument  ensure  air  tight  seals.  The  end  caps  are 
fitted  inside  5.5  in.  pipe  to  form  support  for  the  electronics,  and  to  form  a  flat  base  for  the 
bottom.  A  small  ring  around  the  outside  in  the  middle  of  the  tube  is  used  to  support  the 
device  in  an  arbor  press  for  assembly  and  disassembly.  Two  test  ports  to  in  the  fore  and 
backing  volumes  are  closed  with  short  lengths  of  stoppered  Tygon  tubing  when  not  used 
for  testing.  A  divider  in  the  center,  machined  from  4  in.  PVC  rods  stock,  separates  the  fore 
and  backing  volumes. 

In  order  to  keep  the  sensor  at  a  constant  ambient  temperature  and  maintain  short 
plumbing  connections  (for  good  high  frequency  response),  the  diaphragm  assembly  was 
separated  from  its  electronics  and  mounted  inside  the  fore  volume  on  a  short  length  of  brass 
tubing.  To  prevent  dirt  and  moisture  from  accumulating  in  the  sensor  the  fore  volume  side 
of  the  sensor  is  protected  with  a  short  right  angle  connector.  Wires  from  the  sensor 
penetrate  the  top  cap  through  a  1/8  in.  pipe  fitting  filled  with  a  sealing  compound.  Once  the 
sensor  is  in  place,  the  top  is  filled  with  5  stainless  steel  scrubbing  pads  to  increase  the 
thermal  capacity,  reduce  internal  temperature  gradients,  and  reduce  internal  resonances. 

The  backing  volume  contains  only  the  capillary  tube.  Roughly  27  in.  of  0.03  in. 
capillary  tube  slips  inside  a  short  piece  of  1/16  in.  Tygon  so  that  the  capillary  may  be 
connected  to  the  internal  fore  volume  port.  Once  the  capillary  is  in  place  the  backing 
volume  is  also  loaded  with  5  stainless  steel  scrub  pads. 
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The  bottom  end  cap  is  fitted  with  a  standard  male  garden  hose  connector.  This  is 
left  open  as  the  device  is  assembled  to  prevent  the  backing  volume  from  being 
overpressured  and  damaging  the  diaphragm.  It  is  then  closed  with  a  solid  cap. 

The  top  end  cap  is  more  complex.  It  contains  the  wire  feed  through  for  the 
diaphragm  (4  wires),  a  flexible  hose  to  connect  to  the  input  garden  hose,  a  plumbing 
connector  for  calibration  input  to  the  fore  volume,  the  impulse  pump  calibrator,  and  the 
P305  electronics.  The  impulse  pump  is  normally  connected  to  fore  volume,  but  may 
alternately  be  connected  to  any  of  the  four  connectors  along  the  side  of  the  instrument  to 
allow  impulse  testing  of  the  backing  volume  only  or  combined  testing  of  the  fore  and 
backing  volume. 

Such  a  system  is  expected  to  operate  in  boreholes  at  least  2  m  below  the  surface. 
During  tests  in  the  lab,  the  system  was  placed  in  a  Styrofoam  box  internally  lined  with  thin 
aluminum  plates  and  then  filled  with  Styrofoam  packing  peanuts  (Figure  8). 

During  the  initial  qualitative  tests,  it  was  found  that  the  microbarograph  was 
sensitive  to  small  vibrations  as  a  result  of  flexing  in  the  connector  hose  at  the  top  of  the 
assembly.  Running  the  system  on  a  pier  eliminated  visible  seismic  response.  When  the 
microbarograph  was  closed  and  the  brass  hose  connection  exposed  to  the  ambient  air,  it 
behaved  as  a  gas  thermometer,  registering  the  small  changes  in  temperature  due  to  the 
cycling  building  air  conditioning  This  was  eliminated  by  insulating  the  fittings  with  several 
layers  of  bubble  wrap. 

Calibrator  Design 

The  calibrator  design  had  four  goals.  First,  the  calibrator  had  to  be  remotely 
activated,  ideally  by  using  simple  modifications  to  the  existing  calibrator  electronics  in  the 
digitizer.  Second,  it  must  produce  a  signal  of  sufficient  bandwidth  to  demonstrate  that  the 
dynamic  response  of  the  system  is  unchanged.  Third,  it  should  be  able  to  vary  the 
amplitude  for  varying  noise  conditions.  Fourth,  it  should  calibrate  as  much  of  the  system 
as  practical.  In  this  regard  a  diaphragm  activator  is  less  desirable  than  a  calibrator  that 
produces  a  pulse  in  the  backing  volume.  A  calibrator  which  produces  a  signal  in  the  hoses 
is  even  better.  A  calibrator  which  could  produce  a  known  plane  wave  in  ambient  air  would 
be  ideal. 

Three  different  calibrators  were  examined:  a  modified  bellows  calibrator  from  an 
old  NBS  microbarograph  system,  a  metering  pump  used  as  a  step  generator,  and  a 
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headphone  element  from  a  Koss  high-velocity  open  style  headphone.  Each  of  these  has 
specific  advantages  and  ranges  where  they  proved  useful. 

Bellows  calibrator 

The  bellows  calibrator  (Figure  9)  produces  a  small  volume  change  from  a  DC  gear 
driven  motor.  The  original  calibrator  produced  a  1/3  Hz  signal  when  driven  from  12  Volts, 
close  to  the  high  frequency  end  of  the  original  NBS  microbarograph.  The  motor  was 
replaced  with  a  1  RPM  unit  to  produce  a  constant  displacement  calibrator  that  could  operate 
near  the  desired  low  frequency  comer  of  the  microbarograph.  The  calibrator  is  capable  of 
operating  at  periods  down  to  5  minutes  at  reduced  voltages.  Even  longer  periods  are 
possible  if  the  motor  is  operated  in  pulsed  mode.  Because  the  fore  and  backing  volume  is 
only  1  liter,  the  pump  must  either  be  connected  with  short  tubing,  or  the  volume  of  the 
connecting  tubing  must  be  included  in  the  calculation. 

Impulse  calibrator 

A  surplus  impulse  metering  pump  was  configured  as  a  step  function  calibrator.  On 
each  24  Volt  impulse,  a  piston  moves  a  specific  amount  (adjustable  over  a  small  range 
during  installation),  producing  a  known  change  in  volume.  While  designed  to  work  with 
liquids,  a  film  of  light  oil  allowed  the  pump  to  be  used  with  air.  The  impulse  rise  time  is 
fast,  in  excess  of  0.05  seconds  as  observed  by  the  Gibb’s  effect  in  the  impulse  data.  The 
displacement  may  be  adjusted  from  roughly  0.025  cm3  to  0.250  cm3.  No  leaks  were 
observed  even  after  repeated  cycling,  allowing  the  pump  to  be  permanently  connected  and 
plumbed.  The  pump  may  be  cycled  at  4  Hz  to  produce  a  large  pressure  pulse  on  systems 
which  have  a  long  period  response  of  many  seconds. 

There  are  several  disadvantages.  The  pump  is  surplus  and  exact  off-the-shelf 
replacements  may  not  be  readily  available.  It  requires  a  24  Volt  supply  to  work  reliably  at 
all  stroke  volumes  requiring  either  a  separate  power  supply  or  limiting  dispalcements  to 
only  the  smaller  volume.  During  trigger-  release  cycles  the  coil  produces  a  strong  back 
EMF  that  couples  into  the  signal  line  as  a  spike.  Finally  when  operated  at  24  Volts,  the 
piston  produces  a  strong  vibration  that  becomes  part  of  the  diaphragm  signal. 

Headphone  element 

A  high  velocity  headphone  element  from  a  Koss  open  style  headphone  was  used  as 
the  driver  for  a  calibrator  built  from  plastic  plumbing  (Figure  10).  The  front  of  the 
headphone  element  was  cemented  to  a  reducing  string  and  caulked  to  assure  an  air  tight  fit. 
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The  backing  volume  was  packed  with  foam  rubber.  A  small  compensation  bleed  tube 
connects  the  backing  volume  to  the  fore  volume  through  a  calibrated  orifice  (0.003  in.). 
During  use,  the  element  is  driven  from  a  signal  generator  or  noise  generator  to  produced 
high  to  mid  frequency  drive  signals  (1/10  Hz  to  200  Hz). 

The  small  closed  backing  volume  produces  nearly  constant  displacement  with 
voltage  signal,  since  the  spring  compliance  of  the  headphone  diaphragm  is  dominated  by 
the  compliance  of  the  air  in  the  backing  volume.  The  small  backing  volume  results  in  a 
resonant  frequency  above  500  Hz  in  contrast  to  the  less  than  100  Hz  free  air  resonance. 
Coupling  die  system  to  a  1  liter  input  does  not  change  the  resonant  frequency  significantly, 
indicating  that  the  impedance  of  the  system  is  controlled  by  the  calibrator  backing  volume 
and  bleed  system  rather  than  the  input  impedance  of  the  microbarograph.  Of  course  driving 
a  lower  impedance  system  would  result  in  loading  of  the  calibrator. 

Drive  signals  with  periods  as  low  as  100  second  square  waves  were  used  to 
confirm  that  the  calibrator  produced  a  nearly  constant  displacement  output.  Unfortunately 
we  found  that  the  residual  solvent  in  the  PVC  pipe  cement  would  in  time  attack  the 
diaphragm  resulting  in  small  holes  and  perforations.  While  the  system  could  produce 
strong  enough  signals  to  calibrate  the  microbarograph  itself,  it  did  not  produce  sufficient 
drive  to  calibrate  long  sections  of  porous  hose  at  low  frequencies. 

Microbarograph  Noise  Tests 

Microbarograph  noise,  noise  other  than  that  related  to  a  pressure  change  at  the  inlet 
of  the  microbarograph,  is  the  result  of  electrical  noise,  varying  magnetic  fields,  vibration 
and  tilts,  temperature  strains  within  the  sensor,  temperature  strains  within  the 
microbarograph  housing,  convection  cells  within  the  fore  and  backing  volumes,  and 
mechanical  compliances  within  the  system.  In  the  current  design,  electrical  noise  includes 
sensor  electronics  noise  and  two  additional  sources;  the  single  ended  input  which 
introduces  primary  power  supply  noise  into  the  digitizer,  and  insufficient  RFI  protection  in 
the  sensor  power  supply  circuits.  At  some  installations,  the  electrical  noise  dominates  all 
other  noise  sources.  Varying  magnetic  fields  displace  the  diaphragm  and  change  the 
inductance  of  the  sensor  coils.  Although  the  diaphragm  is  reasonably  stiff,  it  will  act  as  an 
accelerometer  and  vibrations  and  tilts  appear  as  small  signals.  Vertical  vibration  sensitivity 
is  minimized  by  orientating  the  plane  of  the  diaphragm  vertically,  but  this  mayim^pc  the 
sensitivity  to  small  tilts  and  rocking. 

Although  the  sensor  is  constructed  of  a  single  material  to  minimize  temperature 
stresses,  temperature  gradients  across  the  sensor  block  due  to  ambient  gradients  and  heat 
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dissipated  by  the  sensor  electronics  also  introduce  stresses.  Fortunately  the  block’s  large 
thermal  mass  makes  this  only  a  significant  nosie  source  at  low  frequencies.  The  largest 
potential  noise  source  is  caused  by  temperature  changes  and  gradients  within  the 

microbarograph  backing  volume.  A  change  of  less  than  0.0003°C  will  produce  an 

equivalent  input  signal  of  0.1  Pa.  Dissipation  in  any  electronics  within  the  backing  or  fore 
volume  increases  this  noise.  Assuming  that  the  thermal  mass  of  the  system  is  500  gm  H20 

equivalent,  it  would  require  about  1  mW  to  raise  the  temperature  0.00003  °C  in  1  minute. 
Since  the  electronics  dissipate  96  mW,  it  is  essential  to  locate  them  outside  the  fore  or 
backing  volumes.. 


Changes  in  the  backing  volume  dimension,  as  the  result  of  creep,  stress  or  vibration 
create  noise  as  the  effective  volume  changes.  The  dimension  change  required  to  produce  a 
0.1  Pa  signal  may  be  calculated  using  the  gas  relation: 


AV  = 


A PV 
A P  +  P 


where  P  and  V  are  the  base  pressure  and  volume,  AP  is  the  change  in  pressure,  and  AV  is 

the  change  in  volume.  A  volume  change  of  10"6  cm3  or  an  elongation  of  1.2xl0'8  cm  ( a 
strain  of  10'9)  assuming  no  change  in  diameter  would  create  a  0.1  Pa  signal  in  a  1  liter 
backing  volume.  Schedule  40  PVC  is  sufficiently  stiff  that  low  amplitude  seismic  signals 
are  not  a  problem  as  verified  by  loading  the  microbarograph  with  a  1  kg  weight  (equivalent 
to  a  0.1  g  signal)  and  observing  the  lack  of  a  resulting  signal..  For  larger  seismic  signals 
volume  changes  will  be  dominated  by  flexing  of  the  more  flexible  horizontal  hoses  and 
pipes  rather  than  the  hard  PVC  material. 

Three  methods  of  measuring  the  sensor  self  noise  were  considered. 


First,  the  sensor  may  be  isolated  from  the  input  pressure  changes.  For  a 
seismometer  this  is  done  by  blocking  the  mass;  for  a  microphone  it  is  done  by  isolation  in 
an  anechoic  chamber;  for  the  microbarograph,  the  backing  volume  capillary  and  fore 
volume  inlets  are  plugged.  In  the  case  of  a  seismometer  this  technique  overestimates  the 
noise  since  it  is  impossible  to  completely  lock  the  mass  and  laboratories  are  notorious  for 
high  seismic  noise.  In  the  case  of  a  microbarograph,  high  background  pressure  noise  is 
not  a  severe  problem.  First,  while  for  seismometers  ambient  seismic  noise  is  considered 
signal,  but  for  a  microbarograph  ambient  background  noise  coupled  through  the  instrument 
is  noise.  In  contrast  to  the  seismometer,  the  microbarograph  itself  is  supposed  to  have  no 
response  to  pressure  changes,  all  signals  are  introduced  through  the  pipe  array.  Second, 
the  dynamic  range  for  seismic  signals  is  higher  than  that  for  microbarograph  signals.  The 
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advantage  of  ah  isolation  technique  is  that  it  requires  only  one  sensor,  is  quick,  and  gives  a 
reasonable  qualitative  result.  To  the  extent  that  laboratory  noise  couples  through  the  closed 
instrument,  the  technique  will  overestimate  self  noise. 

The  second  method  substitutes  a  dummy  input  for  the  sensor.  For  a  passive 
electromagnetic  velocity  seismometer,  a  resistor  is  typically  substituted  for  the  main  cod. 
For  the  Validyne,  a  matched  set  of  resistors,  a  center  tapped  transformer,  or  matched 
inductors  could  be  substituted  for  the  sensor.  Alternately  a  high  pressure  diaphragm  could 
be  substituted  for  the  ultra-low  pressure  diaphragm.  Since  the  technique  estimates  the 
electronics  noise  only,  it  underestimates  the  total  system  noise,  but  verifies  that  the 
electronic  noise  is  below  the  system  noise  requirement  and  that  stray  RFI  or  EMI  noise  is 
not  significant. 

The  final  method  is  the  signal  coherence  method  which  measures  system  noise 
during  signal  conditions.  Its  primary  disadvantage  is  that  it  requires  two  matched  systems. 
In  a  prior  experiment,  analysis  of  sensor  coherencedetermined  that  self  noise  at  1/50  Hz 
was  20  dB  above  that  at  1/10  Hz.  Since  only  one  system  was  available  in  this  study,  it  was 
thus  impossible  to  measure  coherence.  . 

Noise  Tests  Results 

Four  simple  tests  were  made  to  estimate  the  background  noise  inherent  in  the 
sensors  and  associated  electronics  (Figure  11).  The  microbarograph  was  wrapped  in 
bubble  wrap,  placed  in  an  insulated  box  on  an  isolated  pier  and  allowed  to  temperature 
stabilize  for  3  hours  during  the  late  evening.  The  sensor  was  connected  to  a  24-bit  recorder 
with  the  sensor  and  digitizer  power  supplied  from  a  common  12  Volt  100  amp/hour 
battery.  Data  was  radio  telemetered  to  a  workstation  one  floor  above  and  analyzed  in 
Madab. 

The  fust  test  (Figure  1 1 -dummy  load)  substituted  500Q  resistors  for  the  nominal  20 

mh  sensor  coils.  The  two  resistors  were  glued  in  contact,  insulated  by  wrapping  them  in 
Styrofoam  peanuts,  and  connected  in  place  of  the  sensor  coils.  Data  was  collected  at  40  Hz 
for  12  hours.  The  result  produced  a  wideband  RMS  noise  of  0.006  Pa  or  an  equivalent  35 

uV  RMS  noise.  Small  DC  biases  added  to  either  resistor  result  in  bridge  imbalance  and  can 
create  a  signal. 

The  second  test  (Figure  1 1 -closed  no  bleed)  recorded  data  with  both  the  input  and 
bleed  plugged  resulting  in  a  system  of  two  volumes  connected  by  a  diaphram.  At 
frequencies  below  0.03  Hz,  the  noise  level  rises  as  a  result  of  low  frequency  thermal 
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coupling  or  a  pinhole  bleed  around  the  wire  feed  through.  In  later  tests  we  found  problems 
in  the  seals  of  wire  feed-through. 

The  third  test  (Figure  1 1-closed  system)  recorded  a  sealed  input  port  but  connected 
the  backing  and  fore  volumes  with  a  0.003  in.  bleed  orifice  to  produce  a  high  pass  acoustic 
filter  at  1/60  Hz.  If  the  fore  volume  were  open  or  connected  to  a  low  impedance  source  this 
orifice  would  result  in  a  high  pass  of  1/120  Hz. 

The  fourth  test  (Figure  1 1-HF  closed  system)  is  a  repeat  of  the  previous  test,  but 
with  a  320  Hz  sample  rate.  The  spike  at  60  Hz  is  probably  a  result  of  electrical  noise  rather 
than  magnetic  coupling. 

All  of  the  three  tests  show  similar  results;  a  noise  decay  below  0.3  Hz  proportional 
to  1/frequency  and  above  the  dummy  load  electrical  noise.  The  RMS  noise  from  1/60  Hz 
to  20  Hz  is  0.015  Pa  with  most  of  the  noise  concentrated  in  the  low  end  of  the  spectrum. 
This  is  equivalent  to  88  uV  RMS  noise,  far  below  manufactory  specified  5000  uV.  This 
means  that  sensor  noise  level  will  need  to  be  routinely  verified  for  any  new  sensors,  since  a 
manufacturing  change  could  substantially  increase  noise  levels  but  still  meet  the  advertised 
specification. 

Volume  and  thermal  measurements 

In  order  to  measure  the  effective  fore  and  backing  volumes,  all  test  ports  were 
sealed  and  the  bleed  between  the  two  volumes  was  plugged.  A  syringe  was  then  used  to 
inject  2  cc  of  air  into  one  of  the  test  ports  of  the  closed  fore  and  backing  volume.  Each 
injection  was  repeated  three  times  (Figure  12).  The  temperature  gradient  between  the  fore 
and  backing  volumes  changed  during  the  experiment,  resulting  in  a  background  linear 
trend.  The  overshoot  represents  the  shift  from  adiabatic  to  isothermal  response. 
Exponential  decays  (Figure  13)  were  fit  to  the  detrended  pulses.  This  is  the  classical 
experiment  of  Clement  and  Desmores  to  measure  the  ratio  of  C,/Cv,  the  ratio  of  the  heat 
capacity  at  constant  volume  to  the  heat  capacity  at  constant  pressure.  The  calculated  median 
fore  and  backing  volumes  were  1.07  and  1.16  liters  respectively.  Gamma  (Cp/Cv)  for  the 
fore  and  backing  volume  was  only  about  1 .03  and  the  low  cut  comer  frequency  0.48  Hz 
for  the  fore  volume  and  0.29  Hz  for  the  backing  volume.  The  low  gamma,  less  than  5% 
overshoot,  is  due  to  the  small  volume  and  the  stainless  steel  packing.  This  small  gamma 
means  that  the  effect  of  the  adiabatic  to  isothermal  transition  can  be  ignored  in  computing 
instrument  response. 

2  The  overshoot  is  not  relaxed  to  the  rubber  connectors  or  give  in  the  syringe  since  it  could  not  be  produced  by  squeezing  a  fully 
depressed  syringe. 
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Calibrator  Demonstrations 

Each  of  the  three  calibrators  was  demonstrated  by  performing  one  or  more  simple 
measurements.  The  purpose  was  not  to  exhaustively  test  the  calibrators  but  rather  to 
demonstrate  their  effectiveness  and  deficiencies. 

Impulse  Calibrator 

The  impulse  calibrator  was  connected  between  the  fore  and  backing  volume  and 
adjusted  to  produce  a  3  Pa  pulse.  Tests  in  four  configurations  were  made. 

The  first  test,  where  all  ports  and  bleeds  were  plugged,  verified  that  the  connections 
to  the  impulse  pump  produced  no  new  leaks  and  estimated  the  volume  produced  by  each 
impulse  (Figure  14).  Again,  the  trend  produced  by  the  changing  temperature  gradient  had 
to  be  removed.  The  pump  displacement  was  set  to  the  maximum  that  would  operate 
reliably  at  12  Volts  (the  rated  drive  is  24  Volts).  Based  on  the  backing  volume  calculated 
with  the  syringe,  the  impulse  pump  produced  a  volume  change  of  roughly  1.7  mm3  The 
overshoot  associated  with  the  impulse  is  composed  of  three  components:  a  small 
contribution  due  to  the  gamma  of  the  system;  the  seismic  shock  as  the  piston  hits  the  stop; 
and  the  electrical  and  magnetic  coupling  during  the  rise  and  fall  of  the  solenoid  impulse 
current. 

The  following  tests  fit  simple  one  pole  responses  to  the  recorded  impulse  response 
using  Prony’s  method.  For  cases  where  the  fore  volume  was  open  to  the  laboratory, 
multiple  impulses  were  used  to  create  large  enough  signal  that  could  be  reliably  fit  (Figure 
15). 

In  test  two,  the  fore  volume  was  closed.  The  pump  was  connected  between  the 
fore  and  backing  volume.  A  set  of  two  0.007  in.  and  one  0.003  in.  orifices  were 
connected  in  series  as  the  bleed  resistance  between  the  two  volumes.  When  the  fore 
volume  is  closed,  the  reponse  will  be  related  to  the  equivalent  volume  given  by: 

1 

1  1  • 

- i - 

V  V 

V  fore  Y  back 

Since  fore  and  backing  volume  are  about  the  same,  this  should  result  in  a  period  about  half 
of  that  of  a  measurement  that  leave  the  fore  volume  open  (as  was  done  in  tests  three  and 
four).  In  test  three  (open  fore  volume)  the  same  set  of  orifices  were  used  in  the  bleed 
resistence,  but  the  fore  volume  was  unplugged  and  allowed  to  vent  to  laboratory  air 
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through  a  1  foot  long  5/8  in.  garden  hose.  The  final  test  repeated  test  three,  but  replaced 
the  0.003  in.  orifice  with  a  0.03  in.  by  27  in.  Teflon  capillary  tube. 

Normalized  responses  (Figure  16)  were  determined  and  the  mean  and  standard 
deviation  calculated  for  each  of  the  impulses.  On  the  semilog  plots,  theoretical  exponential 
response  should  plot  as  a  straight  line.  Deviations  are  the  result  of  noise,  such  as  opening 
and  closing  doors  or  elevators  running  during  tests. 

For  the  closed  fore  volume  with  the  orifices,  the  mean  period  was  56.6  seconds; 
for  the  open  fore  volume,  134.2  seconds;  for  the  case  with  the  capillary,  57.6  seconds. 

For  the-final  case  with  the  capillary,  two  standard  deviations  span  periods  from  51  to  66 
seconds.  This  dispersion  could  be  reduced  by  increasing  the  pulse  amplitude  or  selecting 
pulses  with  no  interfering  noise. 

Electromagnetic  Calibrator 

An  electromagnetic  calibrator  using  a  headphone  drive  element  in  plastic  plumbing 
(Figure  17)  was  envisioned  as  a  method  to  make  high  frequency  and  broadband 
calibrations,  for  example  to  determine  the  response  of  the  instrument  up  to  120  Hz. 

The  calibrator  resonance  was  determined  by  observing  the  frequency  where  the 
drive  current  was  minimum.  When  operating  in  free  air  with  the  backing  Volume  sealed, 
the  resonance  was  550  Hz.  Attaching  the  calibrator  to  the  input  of  the  microbarograph  did 
not  change  the  resonance  point  appreciably.  The  calibrator  transfer  function  was 
qualitatively  measured  by  comparing  the  drive  signal  to  the  acoustic  signal  recorded  with  a 
laboratory  B&K  model  4165  microphone  sealed  into  the  front  cavity.  This  produced  a 
response  that  was  almost  flat  from  2  Hz-200  Hz.  Frequencies  below  2  Hz,  the  low 
frequency  comer  for  the  microphone,  were  not  examined. 

The  calibrator  was  sealed  to  the  microbarograph  (Figure  17)  and  driven  with  a  14 
second  period  square  wave  at  ±1.5  Volts,  ±0.47  Volts,  and  ±0.15  Volts  peak-to-peak. 

This  was  done  once  with  the  calibrator  backing  volume  open  and  once  closed  (Figure  18). 
With  the  calibrator  backing  open,  there  is  an  obvious  distortion,  particularly  at  the  highest 
drive  level.  This  could  represent  a  diaphragm  distortion,  a  leak  opened  during  high  drive 

levels  or  even  a  thermal  effect  as  the  voice  coil  heats.  At  ±0. 15  Volts  drive,  with  an  open 

backing,  there  is  a  high-cut  effect  as  well.  Signals  with  the  backing  volume  closed  are 
substantially  smaller  and  undistorted.  The  use  of  a  small  sealed  backing  volume  limits  the 
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diaphragm  movement  to  a  linear  zone  and  results  in  a  calibrator  with  a  constant 
displacement  response. 

Using  the  above  results,  the  peak-to-peak  volume  was  calculated  for  each  drive 
level  with  and  without  the  closed  backing  (Figure  19).  The  volume  displacement  is  linear 
with  drive  level  in  both  cases,  but  the  displacement  is  9.7  times  smaller  for  a  closed 
backing  volume.  This  seems  reasonable  based  on  the  ratio  of  the  calibrator  backing  volume 
to  the  microbarograph  fore  volume.  With  a  closed  backing  volume  the  calibrator  produces 
a  displacement  of  about  0.1  cm3,  sufficient  to  calibrate  the  microbarograph  with  a  closed 
fore  volume,  or  to  inject  a  signal  into  the  backing  volume,  but  insufficient  to  drive  a  system 
of  hoses.  The  drive  efficiency  (Figure  20)  shows  the  nonlinearity  exhibited  by  the 
unbacked  system  compared  to  the  backed  system. 

'It  is  impossible  to  operate  the  calibrator  for  extended  times  with  a  sealed  backing 
since  a  small  ambient  pressure  change  will  change  the  equivalent  spring  constant  or  rupture 
the  diaphragm.  To  avoid  this,  a  0.003  in.  orifice  was  used  to  connect  the  headphone 

backing  to  the  drive  output.  A  100  second,  ±0.45  Volt,  square  wave  was  recorded  by  the 

microbarograph  operating  with  a  sealed  backing  volume  (Figure  21).  Prony’s  method  was 
used  to  fit  a  one  pole  RC  filter  to  determine  the  comer  frequency.  With  the  0.003  in. 
orifice  the  calibrator  had  a  low  cut  comer  of  1/91  Hz,  low  enough  to  calibrate  the 
microbarograph  over  the  full  bandwidth. 

To  demonstrate  the  calibrator’s  wideband  functionality,  four  swept  sine  wave 
sequences  were  recorded  on  an  operationally  configured  microbarograph  (capillary  leak 
installed  in  both  the  microbarograph  and  calibrator)  (Figure  22).  Amplitudes  vary 
depending  on  the  spectral  density  of  the  sweep,  a  function  of  the  time  duration,  bandwidth, 
and  amplitude,  but  it  is  obvious  that  the  amplitude  is  reasonably  constant  over  several 
octaves. 

A  high  amplitude  broadband  calibrator  may  be  used  to  monitor  the  state-of-health  of 
a  porous  hose.  One  end  of  the  hose  was  connected  to  the  microbarograph  and  the  other  end 
connected  to  the  calibrator  driven  by  a  step  sweep  from  20  Hz  to  2  Hz  in  10  increments 
(Figure  23).  The  hose  resonances  appear  as  peaks  at  10  and  8  Hz.  The  hose  was  damaged 
with  a  small  cut  halfway  along  its  50  ft.  length.  A  swept  calibration  resulted  in  a  reduced 
and  differently  shaped  response.  Repair  by  taping  the  hole  shut  restored  the  spectra  to  its 
original  shape  and  amplitude. 
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Bellows  Calibrator 

The  bellows  calibrator  was  briefly  tested  to  determine  the  displacement  volume  of 
the  system  and  the  range  of  frequencies  which  it  could  produce  as  well  as  the  amount  of 
harmonic  distortion  in  the  mechanical  drive.  The  bellows  calibrator  was  used  to  determine 
the  response  of  the  previous  microbarograph  design.  In  this  experiment  the  motor  was 
replaced  with  a  1  RPM  version  to  produce  a  signal  close  to  the  instrument  low  frequency 
comer.  Unfortunately  the  motor  was  destroyed  by  an  overvoltage  before  a  complete  series 
of  data  could  be  collected. 
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Illustrations 


Figure  1.  Skematic  of  an  absolute  pressure  microbarograph  (from 
Haak  and  Wilde,  1996). 


Figure  2.  Skematic  of  a  differential  pressure  microbarograph  (from  Haak 
and  Wilde,  1996). 
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Microbarograph  Tests 
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Figure  4.  Assembled  microbarograph. 
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Microbarograph  Tests 


Figure  5.  Open  microbarograph.  Upper  right  cap  is  the  top 
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Microbarograph  Tests 
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Figure  6.  Inside  of  fore  volume.  A  spare  Validyne  P305D  sensor  is  shown  in  the  lower 
left.  The  Validyne  sensor  block  is  visible  inside  the  fore  volume.  The  stainless  steel  pot 

scrubbers  fill  the  inside. 
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Microbarograph  Tests 
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Figure  8.  Microbarograph  in  insulated  cover  used  during  laboratory  testing. 


Figure  10.  Headphone  calibrator 
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Figure  11.  Noise  tests  of  the  P305D  installed  in  microbarograph. 
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2  cc  volume  change  in  fore  volume 


Figure  12.  Timeseries  from  volume  measurements. 
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Microbarograph  Tests 
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Figure  13.  Exponential  decays  fit  to  the  overshoot. 
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Pressure  (ubar) 


Microbarograph  Tests 
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Impulse  Pump  Calibration  Closed  System 


Figure  14.  Impulse  pump  calibrator  volume  determination. 


Ill 


Pressure  (ubar) 


Microbarograph  Tests 


2/4/9 


112 


Delta  Time  (seconds) 

Figure  16.  Exponential  response  curves  to  impulse  calibration..  Fits  to  test  1  as  solid  line, 
to  test  2  as  dashed  line,  and  to  test  3  as  dashed-dot  line. 
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Figure  19.  Displacement  response  of  the  EM  calibrator. 
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Microbarograph  Tests 
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Calibrator  w/  bleed  into  closed  sensor  +/  0.452  volts 
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.  EM  calibrator  demonstration  sweep  signals . 
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Figure  23.  EM  calibrator  into  the  end  of  a  porous  hose. 
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APPENDIX  5.  INFRASOUND  RECORDINGS  OF  BOLIDES  AND 

EXPLOSIONS 


Jessie  Bonner,  Paul  Golden,  and  Eugene  Herrin 

Southern  Methodist  University  operates  a  four-element  acoustic  array  collocated  with  a 
ten-element  seismic  array  (TXAR)  near  Lajitas,  Texas.  During  the  Fall  of  1997,  three 
events  in  West  Texas  and  New  Mexico  created  infrasound  signals  that  were  recorded  at 
TXAR  and  processed  by  analysts  at  SMU.  These  events  include  a  bolide  that  entered 
the  atmosphere  near  El  Paso,  Texas  and  two  large  surface  (High  Energy)  explosions 
from  the  White  Sands  Missile  Range  (WSMR)  in  south-central  New  Mexico.  For  each 
event,  no  seismic  signal  was  recorded  at  TXAR.  This  fact  combined  with  the  size  of  the 
infrasound  signals  and  a  backazimuth  coinciding  with  the  WSMR  could  bring  concerns 
about  the  resumption  of  atmospheric  nuclear  testing.  In  the  following  pages,  the  results 
from  processing  these  three  events  along  with  plots  of  the  infrasound  arrivals  as 
recorded  at  TXAR  are  presented. 


El  Paso  Meteor 


On  Thursday,  October  9,  1997,  a  series  of  sonic  booms  were  heard  and  a  fireball  (Figure 
5-1)  was  observed  by  many  of  the  residents  of  El  Paso,  Texas  and  surrounding  areas. 
The  cause  of  these  events  was  a  meteorite  that  entered  the  atmosphere  over  West  Texas 
and  southern  New  Mexico.  At  least  four  different  inffasound  arrivals  from  this  event 
were  processed  at  TXAR  (Table  I)  including  the  arrivals  shown  in  Figures  5-2  and  5-3. 


Figure  5-1.  Photograph  of  the  El  Paso  Meteor  (Copyright  1997  Bobby  Boyd) 
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Figure  5-3.  Infrasound  arrivals  (19:11:38.650  and  19:11:41.5)  from  the  October  9, 
1997  El  Paso  meteor  as  recorded  on  the  TXAR  acoustic  array.  The  traces  were 
aligned  to  a  slowness  of  2.710  sec/km  and  a  backazimuth  of  323.848°. 
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White  Sands  Explosive  Tests 


Explosives  tests  were  carried  out  on  the  White  Sands  Missile  Range  on  November  12, 
1997  and  November  19,  1997.  Each  test  consisted  of  detonating  C-4  explosives  stacked 
on  the  ground  surface,  thus  coupling  of  the  seismic  energy  was  limited.  This  is  evident 
by  the  fact  that  seismic  energy  did  not  propagate  to  regional  distances,  and  thus  no 
seismic  signals  were  recorded  at  TXAR  for  either  blast.  Infrasonic  signals  were 
recorded  for  each  blast  and  the  results  of  processing  the  data  is  shown  in  Tables  II  and 
m.  The  signals  from  the  November  12  explosion  had  less  signal-to-noise  ratio  than  the 
November  19  test  (as  evident  in  Figures  5-4  through  5-7)  and  were  difficult  to  process. 


White  Sands  Explosion#l  12-Nov-97  10000  lb.  C-4  surface 

Location:  Lat  33.6205  Long  106.4797 
Distance  to  TXAR:  552.2  km 


Backazimuth  from  TXAR:  331.3° 
Origin  time:  17:47:00.0248 _ 


Infrasound  Arrivals 
(UTC) 

Backazimuth 

Phase  Velocity 
(m/sec) 

U  (m/sec) 

18:16:27.850 

332.5 

338 

313 

18:17:46.150 

335.3 

341 

299 

18:18:29.050 

336.9 

342 

292 

Table  EL  Infrasound  arrivals  from  the  White  Sands  test  on  November  12, 1997. 


White  Sands  Explosion#2  19-Nov-97 
Location:  Lat  33.6209  Long  106.4797 
Distance  to  TXAR:  552.2  km 
Backazimuth  from  TXAR:  331.3° 

Origin  time:  18:00:00 _ 

Infrasound  Arrivals  Backazimuth 

_ (UTC)  _ 

18:25:59.375  329.7 

18:28:41.1  332,5 

18:30:29.475  334.2 


20000  lb.  C-4  surface 


Phase  Velocity 
(m/sec) 

U  (m/sec) 

340 

341 

346 

321 

351 

302 

Table  III.  Infrasound  arrivals  from  the  White  Sands  test  on  November  19, 1997. 
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Time  (seconds)  From  18:18:22  UTC  on  November  12, 1997 


Figure  5-4.  Infrasound  arrival  (18:18:29.050)  from  the  November  12,  1997  White 
Sands  explosives  test.  The  traces  were  aligned  to  a  slowness  of  2.934  sec/km  and  a 
backazimuth  of  336.9°. 
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Time  (seconds)  From  18:26:30  UTC  on  November  19, 1997 


Figure  5-5.  Infrasound  arrivals  (times  shown  by  vertical  lines)  from  the  November 
19, 1997  White  Sands  explosives  test  The  traces  are  aligned  to  a  slowness  of  2.941 
sec/km  and  a  backazimuth  of  329.65°. 
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Acoustic 


Time  (seconds)  From  18:26:30  UTC  on  November  19, 1997 

Figure  5-6.  Infrasound  arrival  (18:26:59.375)  from  the  November  19, 1997  White 
Sands  explosives  test  The  traces  are  aligned  to  a  slowness  of  2.941  sec/km  and  a 
backazimuth  of  329.65°. 
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Dispersive  Infrasound  Arrival 


Prior  to  the  November  19,  1997  White  Sands  test,  characteristics  for  infrasound  signals 
recorded  at  TXAR  fell  into  two  categories:  sharp,  impulsive  types  and  “ringy”  signals 
dominated  by  one  or  two  frequencies.  The  18:26:59.359  arrival  from  the  20000  lb. 
White  Sands  test  (Figure  5-6)  is  the  first  noted  infrasound  signal  at  TXAR  exhibiting 
dispersion.  The  signal  is  normally  dispersed  between  periods  of  0.4  and  5  seconds 
(Figure  5-7).  Since  the  start  time  for  this  event  is  known,  a  group  velocity  dispersion 
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Figure  5-7.  Processing  results  for  the  dispersed  infrasound  arrival  from  the 
November  19,  1997  White  Sands  test.  (Upper  Left)  Raw  waveform  showing  the 
dispersed  signal  at  approximately  55  seconds  into  the  record.  (Upper  Right) 
Dispersion  relations  for  phase  and  group  velocity  as  obtained  from  MFA  and  PMF 
analysis.  (Lower  Left)  Extracted  signal  from  PMF.  (Lower  Right)  Frequency 
domain  representations  of  the  raw  data  and  the  PMF  extracted  signal. 


curve  was  established  using  Multiple  Filter  Analysis  (MFA).  This  dispersion  curve 
served  as  the  input  to  the  Phase  Match  Filter  (PMF)  technique  which  allowed  the 


126 


extraction  of  the  dispersed  acoustic  signal  in  the  time  and  frequency  domains  (Figure  5- 
7).  The  cause  for  the  dispersion  is  under  investigation. 
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